FACILITY FORM 602 



r 




L 


J 


Bellcomm 



I Reproduced by 

NATIONAL TECHNICAL 
! INFORMATION SERVICE 

1 Springfield, Va 22151 


Bellcomm, Inc. 

Washington , D . C . 


TR-70— 


6 * 0-1 


SOLAR POINTING VARIATIONS 
IN EARTH ORBIT AND 
THE .IMPACT ON MISSION DESIGN 

February 11, 1970 


Prepared by: 
B. D. Elrod 


Work performed for Manned 
and Space Administration, 


Space Flight, National Aeronautics 
under Contract NASW-417. 



BELLCOMM, INC. 


TABLE OF CONTENTS 


ABSTRACT 


Page 

iv 


1.0 INTRODUCTION 1 

1.1 Solar Pointing Angles (8,t|>) 2 

1.2 System Performance Factors 3 

2.0 DEFINITIONS AND NOTATION 5 

2.1 Matrix Transformations 5 

2.2 Definition of Sun Line and Inertial Reference System ... 6 

2.3 Orientation of Orbital Reference System 8 

2.4 Earth Oblateness Effect ‘. 9 

3.0 DERIVATION OF SOLAR POINTING ANGLES 11 

4.0 GENERAL CHARACTERISTICS OF 8 AND i|) 13 


4.1 

4.2 

4.3 



17 

19 

23 


4.4 

4.5 


4. 3a ijj. 

C 

4 . 3b ip 


$ Variation in Sun-Synchronous Orbits 
Summary of 8,ip Properties 


24 

27 

29 

31 


5.0 EVALUATION OF SYSTEM PERFORMANCE FACTORS FOR 

CIRCULAR ORBITS 33 


5.1 Gravity Gradient Torque- Impulse in Solar-Inertial Mode . * 34 


5.2 Orbital Sunlight Interval 37 

•5.3 Electrical Energy from Solar Arrays 40 

5.3.1 Evaluation of Electrical Energy . 41 

5.3.2 Electrical Energy in Local Vertical Modes 45 



BELLCOMM, INC. 


1X1 


Page 

5. 3.2.1 Solar Arrays Operated with 

Spacecraft Roll (^0) 46 

5. 3. 2. 2 Solar Arrays Operated with Fixed 

Spacecraft Roll (^0) 51 

5.3.3 Electrical Energy m Local Horizontal Modes ... 57 

5. 3.3.1 Solar Arrays Operated with 

Spacecraft Roll (^0) 58 

5. 3. 3. 2 Solar Arrays Operated with Fixed 

Spacecraft Roll 62 

5.3.4 Electrical Energy m Perpendicular- to- 

Orbital-Plane Modes 67 

5.3.5 Summary and Comparison of Optimum Solar 

Array Performance 68 

5.4 Incident Thermal Energy on Spacecraft 73 

5.4.1 Cylindrical Surface Configurations in 

Earth-Oriented Modes 75 

5.4.2 Flat Surface Configuration m a Solar- 

Oriented Mode 80 

6.0 AVERAGE EFFECTS OF THE 3 VARIATION 82 

6.1 Mission Averages 82 

6.2 Mission Average of System Performance Factors 86 

6 . 3 Approximate Evaluation of Mission Averages of 

System Performance Factors 95 

7.0 SUMMARY AND CONCLUSIONS 99 


APPENDIX A Location of Sun Line and Ascending Node at Launch . . . A-l 

APPENDIX B Gravity Gradient Torque m Solar-Inertial Mode .... B-l 

APPENDIX C Evaluation of Electrical Energy from a Solar Array 

on Spacecraft m LV, LH and POP Attitude Modes .... C-l 

APPENDIX D Pointing Angles for Sun-Oriented Solar Arrays and 

Experiments on Spacecraft in LV, LH and POP Modes . . . D-l 

APPENDIX E Sunlight Interval in Elliptical Orbits E-l 


APPENDIX F List of Symbols 


F— 1 


APPENDIX G List of Tables and Figures 


G-l 


REFERENCES 



'Bellcomm, Inc. 


ABSTRACT 


The solar vector orientation relative to an earth 
orbit varies with time due to the earth's motion about the 
sun and motion of the orbital plane due to earth oblateness. 

In this report a detailed description of this variation is 
given in terms of two pointing angles: 3, the minimum angle 
between the solar vector and the orbital plane and \p r the 
location of orbital noon relative to the ascending node. 

Typical time variations are shown and various properties of 
the pointing angles are expressed as a function of orbital 
parameters for both posigrade and retrograde orbits. Specific 
results are included for 35° and 50° , 230 NM circular orbits, 
which have been considered for AAP Workshop missions. 

The solar pointing angles help define important 
spacecraft systems requirements, such as experiment gimbaling 
and solar array articulation in orbit-oriented modes and 
spacecraft maneuvers in a solar-inertial mode. They are also 
significant from a mission design standpoint, since the solar 
vector orientation is related to factors affecting spacecraft 
system performance. This includes, for example, the incident 
thermal energy per orbit on spacecraft and electrical energy 
output from fixed solar arrays on spacecraft in orbit-oriented 
modes. Although related to both 3 and i|j in arbitrary elliptical 
orbits , performance factors are usually a function only of 3 in 
circular and slightly elliptical orbits. A number of factors 
are considered in this work and the relationship with 3 is 
determined. Spacecraft models and attitude modes are chosen 
to illustrate typical properties and effects of the 3 variation, 
as they may affect AAP and future space station missions. All 
results are given in normalized form to facilitate application 
to particular spacecraft configurations. 

System designs are often based on worst case 3 
conditions. Since these occur only at discrete times in a 
mission, such designs may impose excessive requirements or 
may not utilize available performance margins. Consequently 
mission averages of system performance factors are relevant 
to mission planning and are developed in this report as a 
means of assessing the overall impact of the 3 variation and 
arbitrary launch time and date. 
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SOLAR POINTING VARIATIONS IN EARTH 
ORBIT AND THE IMPACT ON MISSION DESIGN 


1.0 INTRODUCTION 


Solar pointing is a prime concern in manned space- 
flight mission planning because of interrelated requirements 
for earth and solar viewing experiments, electrical power from 
solar arrays, and spacecraft thermal control. Specific point- 
ing requirements are governed by the choice of flight attitude. 
For AAP and future space station missions flight attitudes have 

(1-4) 

been considered' for holding: 

(1) spacecraft including experiments and solar 
arrays in a "solar-inertial" mode,* or 

(2) experiments and solar arrays in a solar 

orientation and the spacecraft in an orbit- 
oriented mode such as: LV, LH or POP.** 

Of interest in these modes are the solar pointing 
angles which describe the orientation of the solar vector 
relative to an earth orbit. These vary with time because of the 
earth’s motion about the sun and continuous motion of the orbi- 
tal plane due to earth oblateness. The pointing angles are 
significant since they, are related to certain spacecraft system 
requirements, such as experiment gimbaling and array articula- 
tion in orbit-oriented modes and spacecraft maneuvers in a solar- 
inertial mode. The pointing angles are also useful from a mis- 
sion design standpoint, since the solar vector orientation is 
related to various system performance aspects (e.g., the incident 
thermal energy per orbit on spacecraft and the degradation in 
electrical energy output from fixed solar arrays in LV, LH and 
POP modes) . 


*The "solar-inertial" mode is depicted in Fig. (1-1) for the 
AAP cluster configuration (Saturn-V Workshop; ATM/CSM/OWS) . As 

in the original concept ^ , the ATM telescope is pointed at the 
sun with the CSM/OWS roll axis maintained essentially in the 
orbital plane and normal to the sun line. 

**The attitude modes LV, LH and POP are defined by the orien- 
tation of the spacecraft roll axis which is : parallel to the 

Local Vertical in LV, in the orbital plane and parallel to the 
Local Horizontal in LH, and Perpendicular-to-the-Orbital-Plane in 
POP. 
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The objective of this report is to describe the 
behavior of solar pointing variations and to examine the effect 
on pointing requirements as well as the impact on various sys- 
tem performance factors affecting mission design. In the 
remainder of this section the pointing angles are defined and 
system performance factors considered in subsequent work. are 
stated. 



Figure (1-1) - Solar Pointing Geometry in Solar-Inertial Mode 


1. 1 Solar Pointing Angles (8,iJ>) 


Pointing in an arbitrary direction is usually 
achieved by two rotations (pointing angles) about two ortho- 
gonal axes. In this report two pointing angles (6,if>)* are 
defined relative to the orbital plane, specifically the geocen- 
tric coordinate system ( x n 'Y n ' z n ) shown in Fig. (1-1) . In this 

coordinate system, y n is normal to the orbital plane, z n is 
directed to the orbit ascending node and (not shown) com- 
pletes the right hand system. The angle ip is a rotation in the 


*The angles $ and ip are synonymous with the pointing 

angles n and n used in Reference 5. 
x y 
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orbital plane about y n and & is a rotation out of the orbital 

plane about the axis x^, which lies in the orbital plane, but 

normal to the sun line. Geometrically, 6 represents the angle 
between the orbital plane and the sun line as measured in the 
noon meridian plane*, whereas i|> represents the angle between 
the ascending node (z^) and the noon meridian plane as measured 

in the orbital plane. 

Due to the earth's motion about the sun and motion of 
the orbit ascending node (z ) in the earth's equatorial plane, 

^ . -t 

3 and vary with time. In subsequent work this variation is 
examined in detail and various properties- of 6 and ^ are evalu- 
ated as a function of orbital parameters. Properties such as 
Islmax, l®linax' +max and Kve are direotl y related to space- 
craft maneuver requirements m the solar-inertial mode and to 
experiment gimbaling and array articulation requirements in 
non- sun- oriented modes. 

1.2 System Performance Factors 

The solar vector orientation relative to the orbital 
plane has an effect on the design of many spacecraft systems 
(attitude control, electrical power, and thermal control) as 
well as many earth and solar viewing experiments. Specific 
factors affecting system performance can be related to the 
solar pointing angles. Significant performance factors include 

(1) Bias gravity gradient torque which contributes 
to angular momentum dumping requirements for 
spacecraft controlled by CMGs in a solar- 
inertial mode, 

(2) Duration of an orbital day which affects 
experiment operations and solar viewing time, 

(3) Incident solar energy per orbit which affects 
the electrical energy available from solar 
arrays , and 

(4) Incident thermal energy per orbit which 
affects spacecraft thermal control require- 
ments . 


*The definition of noon is discussed in Appendix A in 
connection with Fig. (A-l) . < 
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For elliptical orbits these particular factors are a function 
of both g and i|>, although has negligible effect for slightly 
elliptical orbits* and no effect for circular orbits. In 
subsequent sections the factors, (1-4) above, are examined in 
detail for circular orbits and a functional relationship with 
3 is determined. In addition, mission average effects of the 
S variation and arbitrary mission launch time and date are 
studied with a view toward assessing certain system requirements 
and performance margins, as discussed in Section 6. Results are 
plotted as functions of mission duration and orbital inclination. 


*A slightly elliptical orbit is interpreted here as having 
an eccentricity <0.01(e.g., a 260 x 200 NM orbit, which is 
within launch capability for the Saturn-V Workshop mission) . 
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2.0 DEFINITIONS AND NOTATION 


In the following subsections mathematical and geometri- 
cal definitions and notational conventions used in deriving 
pointing angles and in subsequent analysis are given. 

2.1 Matrix Transformations 

are related by transform- 


= T Hj (2-1) 


T X =T^) * and u. and u . 

3 

represent a vector expressed in the i and j systems.** The 
matrix T is an elementary transformation if two coordinate 
systems differ by a single Euler angle rotation about a common 
axis . 

The three possible forms are*** 


(2-2) 


(2-3) 


(2-4) 


*The superscripts -1 and f correspond to the matrix inverse 
and transpose operations . 

**Subscripts on a vector denote the system in which it is 
expressed. Frequently abbreviations are used (e.g., u for a 
vector m spacecraft coordinates) . sc 

***Superscripts on transformations indicate axis of rotation 
and subscripts denote the corresponding angles of rotation. For 
brevity the notation eg and sg is always used for the 
trigonometric operations cos CO and sin(g). 


1 

0 

0 

0 

eg 

sg 

0 

-sg 

eg 

eg 

0 

-Sg 

0 

1 

0 

s g 

0 

c g 

c? 

sg 

0 

“S g 

eg 

0 

0 

0 

1 


Different coordinate systems 
ations of the form 



where T is an orthogonal matrix (i.e., 
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The matrix T may be formed from the product of any number of 
elementary transformations corresponding to Euler angles 
y, £,.... with the property that 

T + = (. . . tJ T k T A ) * 

5 v y' 

= (tJ) + (T^) + (T^) f ... = T_J T_^ T_jj... ( j ,k, £=x,y or z) 

(2-5) 


2 . 2 Definition of Sun Line and Inertial Reference System 

To account for the earth's annual rotation about the 
sun it is convenient to define the sun line with respect to the 

geocentric inertial reference system (X ,Y , Z ), henceforth 

Y Y Y 

called the equinox system shown in Fig. (2-1) . Here Y^ is posi- 
tive toward the ecliptic north pole, Z^ is the direction toward 
the sun at autumnal equinox and X completes the right-hand 
system. 


Earth Location 



Figure (2 - 1) - Location of Sun Line Relative to Autumnal Equinox 

+ 

The sun line is defined as the vector S =(0,0,1) m 

s 

the (X , Y , Z ) system shown in Fig. (2-1) where Y is parallel to 
s s s s 

Y and Z is directed toward the sun. The axes Z and Z are 
Y s s y 

displaced by the angle y which is related to the time of year 
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with. y=0° corresponding to the autumnal equinox. In the equinox 
system the sun line is given by 



(2-6) 


Although the earth's orbit about the sun is slightly elliptical, 
the angle y is evaluated here on the basis of a circular orbit 
approximation , * 


Y ~ Y o + Yt 


(2-7) 


where ** 


Y 


360/365.26 = 


0.9856 deg/day 


(2-8) 


is the mean sun line rotation rate, t is current time and y Q is 
the sun line location at t=0. 


*The actual sun line location (y ) can be determined from 

8L 

the orbit angular velocity equation 

Y a = * (1 “ e s )3/2[1 + e s c(Y a~ Y ap )]2 

where e =0.016 is the earth's orbital eccentricity and y =282.38° 
s (12) a P 

is the earth's perihelion' 7 (relative to autumnal equinox). To 

a first approximation in e it can be shown that 

W 


Ay 


e 




iY e = Y-Y a = 2s s U - s( Ya-Y a p n 


Thus, 0i | Ay_ | <0.03 deg/day and 0^(Ay ) <3.8 deg with the 

s max 0 max 

maximum occurring a few days after vernal equinox (when y =192.38°) . 

a 

For the purposes of this analysis the circular orbit approximation 
is sufficient. 

**The earth's period about the sun is 365.2563835 mean solar 
days. 
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2.3 Orientation of Orbital Reference System 

The relative orientation of the orbital and equatorial 
planes and their relationship to the ecliptic plane and sun line 
is shown in Fig. (2-2) for a particular location of the orbit 
ascending node (z n ) . As previously discussed, the nodal system 

(x n ,y n ,z n ) is considered as the reference system for determining 

the solar pointing angles. The Euler angles e,Ji r and i define 

the orientation of the nodal system relative to the equinox 
system (X^,Y^,Z^). The angle e is fixed at 23.45° and defines 

the inclination of the equatorial plane relative to the ecliptic 
plane. The angle defines the location of the orbit ascending 

node relative to Z^. The angle i represents the inclination 

of the orbital plane relative to the equatorial plane. 

As given in Eq. (A-13) of Appendix A, i can be deter- 
mined from 


ci = sa T c£ 0<i<180° (2-9) 

-u “ *“■ 

where a is the launch azimuth measured positively from geograph- 

J-l 

ical north and % is the latitude of the launch site. Evaluation 
of is discussed in the next subsection. , 



Figure (2 - 2) - Relative Orientation of Orbital, Equatorial and Ecliptic Planes 
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2 . 4 Earth. Oblateness Effect 

Two variations in orbital motion produced by earth, 
oblateness are: 1) rotation of the ascending node about the 

equatorial plane and 2) rotation of the line-of-apsides in the 
orbital plane. Motion of the ascending node can be described 
by* 

R = Q + h t (.2-10) 

r ro r 

where R r is the rotation rate, t is current time and R ro is the 

ascending node location at t=0. Similarly, motion of the line- 
of-apsides is described by* 

ft = R + Q t (2—11) 

a po a 

where h is the rotation rate and R is the location of perigee 
a po 

relative to the ascending node at t=0 as shown in Pig. (2-3) .** 

North Pole 


Perigee 


Equatorial Plane 


Ascending Node 


Orbital Plane 

Figure (2 - 3) - Description of Orbital Motion Due to Earth Oblateness 


*This model represents the secular variation due to earth 
6 

oblateness . Periodic variations due to third and higher harmonic 
terms in the earth's gravitational potential function are small 
and may be neglected in this work. 

**Although the line-of-apsides is undefined in circular orbits, 

the orbital plane may be viewed as rotating (in-plane) at the 

rate, fl . 
a 
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The rotation rates 


and ft 
r a 


are given by 


ft 


r 


and 


ft 


a 


where 


ft 


rm 


and 


J 


y 


-ft ci 
rm 


(2-12) 


, c 2. 

4-5 s i 


rm 2 


(2-13) 


Jh 1/2r2 = 9 996 CR / a) ^ 3 . 

a 7/2 ( i -e 2)2 9 ’" 6 (1 _ £ 2 )2 day 


(2-14) 


-3 

1.624 x 10 , dimensionless constant in second 

harmonic term of oblate spheroid 
model of earth’s gravitational 
potential function 

4.6843 x lO 1 ^ NM 3 /day 3 , earth's gravitational 
constant 


R = 3443.9 NM, earth radius 
e = orbit eccentricity 
a = orbit semi-major axis (NM) 

= R+H for circular orbits (e=0) 
H = circular orbit altitude (NM) 


Motion of the ascending node is westward (regressive) 
for 0°<i<90° and eastward (progressive) for 90°<i<180°. Apsidal 
motion is such that perigee advances relative to the ascending 
node for 0°<i<63.45° or 116.55°<i<180° and- reverses for 
63. 45°<i<116. 55° . 
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3.0 DERIVATION OF SOLAR POINTING ANGLES 


The pointing axis (P) to be aligned with the sun line 
can be defined relative to the nodal coordinate system (x n ,y n ,z n ) 

by the pointing angles (g,ij>) described in Section 1.1. To derive 
(B,*) it is only necessary to equate elements of P with corre- 
sponding elements of the solar vector S also expressed in nodal 
coordinates . 


To determine P consider the spacecraft coordinate 

system (x ,y , z ) shown- in Fig. (1-1) which is related to the 
sc sc sc 

(x n/ y n / z n ) system by the transformation 


y = T u = T^ T^f y 

-sc p £ n g -n 

4* 

Since P =(0,0,1) , it follows that 


P = 
-n 




s ^ eg 
-sg 


C if) Cg; 


(3-1) 


(3-2) 


To determine the solar vector S consider the trans- 

— n 

formation between the ( x n 'Y n ' z n ) an< ^ ( X Y »Y^,Z^) coordinate 
systems shown in Fig. (2-2) . 


U a = T n M y = T? T * h y 


(3-3) 


From Eqs . (2-6) and (3-1) it follows that § n is given by 


fs i! 


-n 


=1 S. 


T S 
n - Y 


3 / 


sy(ci ce cO r -si se’)-ci cy 
|-sy(si ce cft r +ci se)+si cy sG 
cy cJ2 r +sy sQ r ce 


(3-4) 
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Equating P n and yields the trigonometric functions 
specifying the pointing angles (B,iji).* 


sg = -S 0 = sy {si ce cfl + ci se) - si cy sfi 


(3-5) 


tarn]/ = — = 

S 3 


^ sy (ci ce cS2 r - si se) - ci cy sS2 j 


cy cJ2 r + sy s^ r ce 


(3-6) 


or in alternate form 


s 0 = ci se sy + si 


i [(iip) S (Y-S r ) - fcp) s (Y+!! r )] 


(3-7) 


tan\J> = 


-si se sy + ci 

/l+ce\ , 

i 2 ) s (IT" 

v -(S 06 ) 

\s (y+« r )j 


[W 

|c(y-G r ) +( 

;i^)c(v« r ) 



(3-8) 


As discussed previously, the angles e and i are fixed whereas 
y and var 

be examined. 


y and vary with time. The effect of this variation will now 


*In general, two solution pairs (.0,^) exist, one pair 
corresponding to -90°£B<90° and the other for 90°£B<270°. In 
subsequent analysis the pair corresponding to -90°<j5<90° is 

used which implies that 3 always represents the minimum angle 
between the sun line and the orbital plane. 
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4.0 GENERAL CHARACTERISTICS OF 3 AND 


The time variation of 3 and is affected by the 
initial values of y and Derivations for y Q and as a 

function of launch conditions are given in Appendix A. From 
Eqs. (A- 3) and (A-5) the initial values of y and are* 


Y 0 = Y<3 l H-l) 

8 ro = I15(t L +t N )-180'>] + a H - S! l (4-2) 

where 

d L - D L + 


= launch time measured from autumnal 
equinox in days (non-integer) 

D l = number of calendar days (integer) 
at launch since September 23 

(00 h :00 U.T. **) 

t^ = launch site Universal Time** at 
launch (0<t t <24.0 hours) 

t = launch site Universal Time** at 
de autumnal equinox (0<t <24.0 hours) 

— cl© — 


and 


*The expression for y Q is based on a circular orbit approxi- 
mation to the earth's motion about the sun. See Footnote *, p.7. 
For more accurate results y o can be evaluated from ephermeris 

tables^, once d^ is determined. 

**See Footnote *, Appendix A, p.A-2 for definition of Univer- 
sal Time (U.T.) at launch site. 
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= 


difference between true noon* and 
apparent (clock) noon at Greenwich 
(0<J t N | ^0. 27 hours) 


As illustrated in Fig. (A-l) of Appendix A, is the launch, site 

longitude relative to the orbit ascending node and is the 

longitude of the equatorial noon meridian plane* relative to 
the autumnal equinox. The expressions given in Eqs.(A-9) and 
(A- 15) for obtaining and & L are 

tan = tany Q ce (4-3) 

tan = tana T SSL (4-4) 


where and SL represent the launch azimuth and launch site 
latitude. The angles and lie in the same quadrants as y Q 
and c* L respectively. ' 

Although these definitions for and y are based on 

3 ro ’o 

instantaneous orbit injection above the launch site, they are 

also applicable for arbitrary orbit injection points or for 

previously established orbits. Only appropriate interpretation 

of parameters is necessary, specifically: <* L is the current 

heading relative to geographical north; % is the current local 

latitude; D T ,t T ,t _ and t. T are times (as defined above) corre- 
L L ae N 

sponding to the current local longitude.. 

To illustrate the nature of 8 and iji as a function of 
time, the solutions of Eqs. (3-5) and (3-6) for 8 and i|> are 
plotted in Fig. (4-1) for a specific set of orbit and launch 
conditions.** Thes'e results are based on a noon launch at autumnal 
equinox (y =0° : t=t =12 . 0 , D=t =0) . The orbit and launch 

O J_i a.e Xj jn 


*See discussion in Appendix A regarding the definition of noon. 
In this formulation true noon corresponds to ^+-^=12.0. 

♦♦Expressions for evaluating the initial pointing angles (B 0 r^ 0 ) 

for arbitrary launch conditions are given in Appendix A. 3. 
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Time (days) 


Figure (4 - 1 ) * Typical Variation in Solar Pointing Angles. 
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parameters are based on a 35°, 230 NM circular orbit (i=35°, 
H=230 NM, e=0) which corresponds to & r =-6. 51 deg/day. The 

effect of launch time on 6 and is shown in Fig. (4—2) where 
the curves correspond to four different launch times 
(t L =6, 12, 18,24) with the other parameters the same as those 

used for Fig. (4-1) . 

The results shown in Figs. (4-1) and (4-2) indicate 

that ifi increases monotonically with time, in fact, almost 

linearly. The curves also indicate that 6 varies harmonically 

within an envelope whose maximum magnitude occurs at the two 

solstices (y =90°, 270°) with the value | @ f ={i+e). It is of 
o max 

interest to examine these properties in general. Other charac- 
teristics of 0 and if> are also of interest, namely 

max f uVs 

and li| max » l n the following subsections analytical expressions 
are obtained which relate these various characteristics to y, 
and i. 


4.1 0 Envelope 


The 0 envelope defines the possible extremes of 6 
which may exist at any time of year (y) . At any particular time 
(y fixed) extreme values of 0 can be determiried by differentiat- 
ing s0 in Eq. (3-5) with respect to fi r , since 


3 S0 = 9 S0 . 30 = 30 

3R 90 9S2 P 9Q 

r r r 


Hence, it follows that 


9 s 0 
9fl 

r 


-si(sy sft r ce + cy c& r ) 


0 


(4-5) 


which implies that 

tan£2 r = -1/ (tany ce) (4-6) 

The two possible solutions for £2 r give rise to the upper and 

lower boundaries of the envelope. Using Eq. (4—6) with Eq. (3-5) 
yields the boundary equations in terms of s0. 

s0 - sa^ ci + si ca^ = s(a^+i) = s 0^ 


(4-7) 
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where 


sa^ - sy se 


(4-8) 


The + and - algebraic signs correspond respectively to the upper 
and lower envelope boundaries. The angle -a^ can be interpreted 

physically as the sun declination relative to the equatorial 
plane. 


The time of year (y) corresponding to the maximum 
points on the envelope, namely |b|__ , can be determined from 

ill cl X 


3S 0 


eb _ 


3s8 


eb 


3d. 


3y 






91 se = 0 
Ca d 


The solutions for y obtained from c(a^±i)=0 yield 

A 

sy = + ci/se = + ci/c (90 °~e) = sy 


(4-9) 


(4-10) 


However, these apply only if (90°-e)^ii (90°+e) . 
solutions for y follow from cy=0 which yields 

y = 90°, 270° = y 


The other 

(4-11) 


This case applies for 0<i<(90°-e) and (90°+e) <i<180° . 

The values of |fH max and the values y=y and & r =S2 r * 
yielding J 3 ] are summarized in Table 4-1 for the four appli- 

ITlclX 

cable ranges of orbital inclination. 


TABLE 4-1 — /? Envelope Parameters as a Function of Orbital Inclination 


Inclination** 

Envelope^5^ e 

Parameters 

(A) 

0°< i < (90° - e) 

<B) ! 

(90° -e)<i<90° 

(C) 

90° < i < (90° + e) 

(D) 

(90° + e) < i< 180° 

1^1 max 

i + e 

CD 

O 

o 

90° 

180° - i + e 

A 

7 

90°, 270° 

± sin 1 (ci/se) 

± sin — ^ (ci/se) 

90°, 270° 

A 


'V 

„ 'V 


n r *** 

0° 


180 + S2 r 

180° 


*S2 r is obtained by solving Eq. (4-6) with y=^. 

**e =23.45° 


***0 = tan -1 [-1/ (tan y ce) ] ; -90°$Q *90°. 

3T 
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Plots indicating the form of the 6 envelope for the 
four ranges of inclination (A,B,C,D) in Table 4-1 and two special 
cases (i=0 o , 180° ) are shown in Fig. (.4— 3} . The upper and lower 
boundaries are inverted images of each other about y=0° or 180°. 
For ranges A and D the upper and lower boundaries are separated 
by a constant angle, 2i and 360°— 2i respectively, as may be 
deduced from Eq. (4— 7) . For ranges B and C this is true only for 
the values of y between successive peaks of opposite sign. The 
The minimum boundary separation in ranges B and C is C180°-2e) 
which occurs at y= 90 ° and 270 °. 



'4.2 



Figure (4 - 3) - Effect of Orbital Inclination on 0 Envelope. 


The angle 0 varies harmonically within the 8 envelope 
except for i=0°,180° when the envelope collapses to a median 
curve defined by 

. -i 

6 e = sin Csy se) (i=0°) (4-12) 

8 e = sin -1 C-sy se) (i=180°) 


(4—13) 
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which, is the sun declination profile relative to the equatorial 

plane. Unless i=0°,180° or a particular inclination (i ) 

s s 

corresponding to a sun-synchronous orbit*, (3 crosses the median 
curve between successive intersections with the upper and lower 
boundaries as illustrated in Fig. C4-4) for ranges A and B of 
Table 4-1. 



Launch Time Noon (AE) 



AE - Autumnal Equinox 
W$- Winter Soisl ice 
VE - Vernal Equinox 
SS - Summer Solstice 


Figure (4 - 4) - Typical (3 Variation for Orbital Inclination in Ranges A and B. 


The time rate of change of B can be calculated by 
differentiating Eq. (3-5) which yields 


d(sB) 

dt 


c3 h 



or 


3 = 


K .Q + 
1 r 

cB 


k 2 y 


(4-14] 


(4-15] 


*See Section 4.4. 
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where 

= -siCsy ce sfi r + cy c^ r ) (.4-16) 

K 2 = cyCsi ce c^ r + ci se) + sy s« r si (4-17) 

The 3 variation shown in Fig. C4-4a) suggests that [b| reaches a 
local maximum as 3 crosses the median curve. A formal solution 

for all | 3 [ obtained by differentiating 3 with respect to 

y and in Eq. (4-15) is a tedious exercise. One solution is 

| 3 | = | (-^ r )si + ys (i+e) | = | ^ rm s i ci + ys ( i+e) | _ (4-18) 

and the corresponding y and are 


; =0°; y = 0° , 180° 

0° 

< i < 90° 

(4-19) 

: = 180°; y = 0°,180° 

90° 

< i < 180° 

(4-20) 


In Fig. (4-5) the dashed curve represents |g| plotted versus 

inclination. The solid curve is a numerical solution for Igl ,,, 

1 1 max 
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obtained by scanning all y, Prom the figure it is evident 

that j 3 1 is indeed the upper limit of 1 3 1 when i is in the range 
A or D of Fig. (4-3). For i in the range B or C, [ 8 | is also the 
upper limit of jg| except in the vicinity of the inclination/ i . 

The upper limit on occurs at an inclination in 

A luciX 

range A. Differentiating | 6 j with respect to i leads to the 
equation 

c2i = (-y/^ rm ) c (i+e) (4-21) 

Applicable solutions for i lie in the range 45° <i< (90'°-e) and 
depend on & rm * ** A plot of the corresponding upper limit on 
I ^ I max s k° wn Fig- (4-6) as a function of altitude for circu- 
lar orbits. For 230 NM the absolute upper limit is | 3 [ = 

4.90 deg/day at an inclination of 46.3°. This compares with 
I ^ I max = deg/day for i=35° at 230 NM as shown in Fig. (4-5) . 



Figure {4 - 6) - Upper Limits of |£j max and Corresponding Orbital Inclination 
vs Altitude (Circular Orbits). 


*Q is a function of the orbit semi-major axis and eccentri- 
rm J 

city in elliptical orbits and a function of orbital altitude alone 

in circular orbits. See Eq. (2-14) . 
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4.3 Properties of ^ 

The monotonic variation of tf* for the case shown in 
Fig. {.4-1} indicates that ip remains positive and fluctuates about 
some average rate, ^ . General conditions under which these 

cLV© 

properties hold can be determined from an expression for i|). 
Differentiating tan ip in Eq. C3-6) with respect to time yields 


d ^ (X • 2 , 13 1 3 

^•(tanifi) = ij; sec ip - ^ (4-22) 

S 3 

where S 1 and S 3 are defined by Eq. {3-4} . Substituting for 
S^, S^, S-^, and and using Eg. (2-12) yields 


where 


ip = 


F 2 rm + 6 V 
c 2 S 



F = ci(ci - sy se s£) 

G = (ci ce - c^ r si se) 


(4-23) 

(4-24) 

(4-25) 


The coefficient F is positive for all when the 

orbital inclination i is in ranges A and D defined in Table (4-1) 
and Fig. (4-3)*. In range A the coefficient G is positive for 
all but in range D it is negative for all & r .** Since & rm 

and y are positive, ^ is always positive in range A and will be 
positive for all i in range D provided that is sufficiently 

large. For circular orbits this is the case for orbital 


*Since the extreme value of sysB is s(i+e) in range A and 
-s (i-e) in range D, it follows that F m ^ n =cicec (i+e) >0 in range A 

and F . =cicec (i-e) >0 in range D. 
mm 

**Extreme values of G occur for & r ~ 0° in range A where 

G . =c(i+e)>0 and for fl =180° in range D where G__ =c (i-e) <0 . 
mm r ^ max 
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altitudes less than 690 NM. * Since coefficients F and G may be 
positive or negative for i in ranges B and C, the monotonic 
variation of ^ will not prevail, but may occur at intervals. 


Although <]> is positive for i in ranges A and D it 
fluctuates about an average rate i ave r as observed previously. 

This average rate can be calculated by observing the number of 
revolutions of over a particular time .interval. Revolutions 
of i/) can be referenced to a point where which occurs when 

S^O in Eq. (3-6) and implies that 

tan& r = -l/(tany ce) (4-26) 


the same condition, Eq. (4-5), for g to intersect the g envelope.** 

I 

Let & r , y' and lj^be values of y and i|> at time, t=0. Suppose 

the ratio |fl r |/y is rational; then at some time t=x=2irk/y, and 
y become 


» . > _ 

£2 — £2 + fir = 0 + 2 7tia 

r r r r (4.27) 

' . ' 

y = y + yx = y + 2 irk (m,k integers) 


where (+) signs account for the sign of *** Inspection of 

solutions to Eq. (4-26) reveals that in the same interval x, tji 
makes exactly m+k revolutions so that 


Thus , 


i|)(x) = 2ir (m+k) + 


(4-28) 




ave 


^ (t) — xp o 
x 


2 it (m+k) 
2 irk/ y 



(4-29) 


*The altitude limit for iJj> 0 is a minimum at # i=90°+e. At 
inclinations above 90°+e the altitude range for tj>>0 expands 
considerably (e.g., altitude limits at i=90°+2e and i=180° are 
2206 NM and 3075 NM) . This limit is found from Eq. (2-14) after 
minimizing g in Eq. (4-23) with respect to y and £2^. and evaluating 

& such that g . =0 . 
rm min 

**Hence, revolutions of $ (relative to mr/2, n=l,3,...) are 
synonymous with successive intersections of g with either the 
upper or lower envelope boundary. 

***In what follows, the upper sign in (+) or (+) applies to 
range A and the lower sign to range D. 
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» « • 


But m/k = + ft r /y from Eq. 04-27} so that 


^ = (+ h + y) = 

ave r - 


f (ft ci + y) , Range A 


U riu ei + Y) / Range D 


(4-30) 


Although based on a rational ratio for |fi r |/y f this 

result is approached in general for t sufficiently large. To 
show this , let 

» » 

^r ~ = + ^ re 2nnr 

, . . (4-31) 

Y = y + y t = Y t 2kfr 

at time t, so that 

)|/(t) = 2tt (m+k) + (4—32) 


where -ir< (ft ) <n represent the residue in ft and i/> for 

re e r 

| ft r I / Y not rational. Thus 


^ ( t ) - 1|> 




ave 



* ^ 

X — 

k 2ir 


(4-33) 


Solving for m/k from Eg. (4-31) and substituting into Eq. (4-33) 
yields 


r ave 



(4-34) 


which approaches the result in Eg. (4-30) as k (i.e.. 


t) increases. 


A plot of ^ for a 230 NM circular orbit and orbital 
ave 0 

inclination in ranges A and D is shown in Fig. (4-7). Since i|> 

may be either positive or negative for i in ranges B and C, 

| Travel generally will be less than the result given in Eq. (4-30) . 

Extrapolation of the ^ ave curves for ranges A and D into ranges 

B and C is also shown in Fig. (4-7) for comparison. 


*Hence, successive intersections of 3 with either the upper 
or lower envelope boundary occur at intervals of 2-n/\p , on the 
average. av 



(deg/ day) 
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Figure (4 - 7) - ^ ava and max vs Orbit Inclination (Circular Orbits). 
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The ripple on the curve of ^ vs t shown in Fig. (4-1) 

reflects the fluctuation of ip about . At orbital inclina- 

ave 

tions where 6 may approach 90°, iji can be large at times , in fact 
infinite when 6=90°. These aspects are illustrated in Fig. (4-8) 
where typical 6 and ip variations are plotted versus time for two 
orbital inclinations : one in range B and one in range A near 

(90°-e) .* 


As is apparent from Fig. (4-8) local maxima of occur 
as 8 reaches the 8 envelope. The upper bounds on if> can be 
shown to be 


/ ® rm ci ce + y 


c (i+e) 


ij) = ( 

r max 1 


4- “ 


(+) Range A, (-) Range D 

(4-34) 

Ranges B and C 


The corresponding y and & r for each range of inclination are as 
given in Table (4-1). For comparison with ip these results are 

3.V6 

plotted versus i in Fig. (4-7) for a 230 NM circular orbit. 


*Having is not a serious problem since it merely repre- 

sents an instantaneous condition for satisfying the solar point- 
ing requirement as the orbital plane passes through an orientation 
where it is normal to the sun line (8=90°). As a practical matter, 
the "jump" in can be accomplished in finite time with spacecraft 
maneuvers or experiment gimbaling without significant time inter- 
ruption of a solar pointing mode. 




AE — Autumnal Equinox 
WS - Winter Solstice 
VE - Vernal Equinox 
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4.4 g Variation in Sun- Synchronous Orbits 

7 

A sun- synchronous orbit has been defined as "an earth 
orbit whose plane rotates about the earth's axis just one revolu- 
tion per year and remains in step^with the earth's rotation about 
the sun."* This implies that Q^=y or in view of Eq. (2-12) that 

ci = -y/9, a ci (4-35) 

* rm ss 

Since y and 2 are positive, sun-synchronous orbits are 
necessarily retrograde orbits (9Q°<i ^180°) . The relationship 

5 5 

between orbital inclination and altitude for attaining sun- 
synchronous circular orbits based on Eqs . (2— 14) and (.4—35) is 
plotted in Fig. (4-9) . For a 230 NM circular orbit the required 
inclination is 97.13°. 



Figure {4 - 9) - Orbital Inclination and Altitude Requirements for Sun-Synchronous Circular Orbits . 


The g variation in a sun-synchronous orbit can be 
determined from Eqs. (2-7) , (2-10) and (3-7) where the nodal 
rotation is given by 

®r = fi ro + ^ ± Y o = ^r + Y (4 " 36) 


and 


J2 — 


ro 


y. 


(4-37) 


is a -constant depending upon initial conditions. From Eqs. (3-7) 
and (4-36) it follows that 


sg = + g sy + h. s (2y+Si r ) (4—38) 

*Since the earth's orbit is slightly elliptical, "remains 
in step" is not true exactly, only on the average. However, the 
deviation is small. 
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where 



14-39) 


g b ci se 


£4-40} 



(4-41) 


Hence, in sun-synchronous orbits B varies about a constant, 

S =sin (-sisfi ), with, two dominant frequency components, y and 

S 2T 

9 

2y. Curves of B vs t are plotted in Fig. (4-10) for various 

values of Q . 

r 



Figure (4 - 10) - /3 Profiles in Sun-Synchronous Orbits . 


These results indicate that B does not make successive 
intersections with the B envelope as in non-sun-synchronous 
orbits. Since the constants g and ft in Eg. (4-38) are small, B 
remains in the vicinity of except near {^=+90° where it 

follows the envelope variation.* The results are similar 


*For fi r =+90°, it follows from Eqs. (3— 5) and (4—6) that 
sB = ci(syse) + si [l-s^y (1-ce) ] 

which is nearly equivalent to the B envelope equations. The first 
term in this equation and that in Eq. (4-7) are equal and the .second 
terms are identical through the first two terms of a series 
expansion in e. 
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for 90° < | & r | <18Q ° with. & varying about zero for G r =+18G°. Sun- 

synchronous orbits may be of particular interest for long dura- 
tion missions where a relatively flat 8 profile is of importance 
in certain photographic or* earth viewing experiments. 

4.5 Summary of 8 , Properties 

General properties of the solar pointing angles (8/tJ>) 
have been evaluated as a function of time of year {.sun line 
location, y) , ascending node location C^ r ) , orbital inclination 

(i) and the sun line and nodal rotation rates (y,fi r ). Specific 

properties depend on the range of orbital inclination (A,B,C, or 
D) defined in Table (4-1). In general 8 varies harmonically within 
an envelope which depends on both y and i. Upper limits for the 
envelope are given in Table (4-1). Limits on 1 8 I within the 8 

envelope are shown in Figs. (4-5) and (4-6) and may be evaluated 
from Eq. (4-18) . 

The angle i|> increases monotonically for i in ranges .A 
and D, but not necessarily in ranges C and D where the orbital 
plane may pass through an orientation normal to the sun line 
( )8 i =90° ) . In all cases where | b| approaches 90°, 1 | increases 
sharply with | ip | max =0 % theoretically. A summary of properties of 

(B,</>) for typical AAP Workshop missions is given in Table (4-2). 


TABLE (4-2) Summary of {/?, ip) Properties for Two 230 NM Circular Orbits) 


n 


^a 

1^1 max 

l»l max 

^ave 

^max 

35° 

-6.5° / day 

9.3° /day 

58.45° 

4.57° / day 

7.5° / day 

13.2° /day 

50° 

-5.1° / day 

4.2° /day 

73.45° 

4.86° / day 

6.1° /day 

19.9° / day 


As noted in Section 4.3a, 8 intersects the 8 envelope 
when ii>=nir/2 (n=l,3,...). Although this does not occur uniformly 
in time, successive intersections with either the upper or lower 
envelope boundary occur at intervals of 360°/ij> ave , on the aver- 
age. Hence, the "average period" of the 8 variation is 48 days 
in a 35°, 230 NM orbit and 59 days in a 50°, 230 NM orbit. 
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In sun-synchronous orbits (& r = y) , 3 does not make 
successive intersections with the B envelope, but varies slightly 
about a constant, 3 = -sin is.s(ft -y )] as shown in Fig. (4-10) . 

S jL 3TO O 

Since g can be selected by appropriate choice of launch condi- 

tions (& ro f y q ) , sun-synchronous orbits may be of interest where 

a relatively flat 3 profile is of importance for photographic and 
earth-viewing experiments , electrical power from solar arrays and 
spacecraft thermal control considerations. 

Physically, tp represents the angle between orbital 
noon and the ascending node as shown in Fig. (1-1) . Hence, noon 
may advance relative to the ascending node by as much as 13.2°/ 
day in a 35°, 230 NM orbit and 19.9°/day in a 50°, 230 NM orbit. 
This is of interest in the Saturn-V Workshop (SVWS) mission 
where the ascending node will be located by ground tracking and 
orbital noon by calculating ijj. 


In maintaining the solar-inertial mode spacecraft 
motion can be described m terms of azimuth and elevation maneu- 
vers relative to the nodal coordinate system (x ,y ,z ) in 

n u n n 

Fig. (1-1) . The elevation maneuver rate (B^) is simply B, but 
because of orbital plane motion (& )*, the azimuth maneuver rate 
(* a ) is 

~ ^ “ ^a (4-42) 


Hence, the maximum azimuth maneuver rate, (i _)__„, is 3.9°/day 

u. max 

in a 35°, 230 NM orbit and 15.7°/day in a 50°, 230 NM orbit. 

The maximum required angular momentum of the SVWS (including 

the effect of o with 3 and i ) can be shown to be** ^3.0 ft lb 
“r a Y a 

sec in the 35° orbit and 'v-lO.O ft lb sec in the 50° orbit, which 
is well within the angular momentum capacity of the CMG system. 


Pointing angles for maintaining solar arrays and experi- 
ments in a solar orientation with spacecraft in orbit-oriented 
modes are related to B. This relationship is developed in Appendix 
D for spacecraft in LV, LH and POP modes. 


*See Footnote **, p. 9. 


* Calculations are based on mass data in Reference 8. 
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5.0 EVALUATION OP SYSTEM PERFORMANCE FACTORS FOR CIRCULAR ORBITS* 


Earth orbital mission design is influenced by a 
variety of conflicting constraints, some due to mission objec- 
tives and others due to orbit environmental aspects. Of 
importance in the latter category, for example, are the effects 
on attitude control, electrical power and thermal control system 
performance requirements imposed by motion of the orbital plane 
due to earth oblateness and motion of the earth about the sun. 
These effects on system performance can be related to the solar 
pointing angles, $ and ifi, although in circular orbits only 3 may 
be involved. 

The objective of this section is to examine various 
system performance factors which are related to 6. Specifically, 
this includes evaluation of: 

(1) the bias gravity gradient torque-impulse on 
spacecraft in a solar-inertial mode, 

(2) the orbital sunlight-interval, 

(3) the electrical energy per orbit available from 
solar arrays in LV, LH and POP attitude modes** 
with incrementally articulated solar arrays , 
and 

(4) The incident thermal energy per orbit on space- 
craft in LV, LH and POP modes** and on a flat 
solar panel in sun-oriented modes. 

In this work spacecraft and solar array models are 
chosen primarily to illustrate general properties and the 
influence of the 3 variation. All results are given in norma- 
lized form which permits simple application or approximation to 
particular spacecraft configurations or changes in orbital 
altitude . 


*Results in this section for circular orbits are essentially 
the same for slightly elliptical orbits of equal period. (See 
Footnote *, p. 4) 

**See Footnote **, p. 1. 
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5. 1 Gravity Gradient Torque- Impulse in Solar-Inertial Mode 

An asymmetric* spacecraft in orbit experiences a 
torque due to the earth's gravitational field unless one princi- 
pal axis remains parallel to the local vertical. In the solar 
inertial mode the orientation of the spacecraft geometric axes 
relative to the sun line and the orbital plane is 

SC sc sc 

defined by the solar pointing angles as illustrated previ- 

ously in Fig. (1-1) . For the Saturn-V Workshop (SVWS) model shown 
Fig. (5-1) the spacecraft principal axes are assumed to be dis- 
placed from the geometrical axes by a single rotation ($) about 

the roll axis (x ).** 

s c 



y sc 



Saturn V Workshop Configuration 
(ATM Solar Arrays Omitted) 


Sun 


Figure {5 - 1) - Spacecraft Orientation in Solar - Inertial Mode . 


*Here asymmetry implies unequal principal moments of inertia 

(i.e. I 7*1 ). 

v xp r yp r zp 7 

**Coincidence of the principal axis of minimum moment of 
inertia (X ) and the roll axis (x ) is a reasonable assumption 

P S C g 

for this analysis. In a current dry-workshop configuration the 
actual displacement is only 3.6°. However/ the rotation of 
transverse principal axes is 23°. 
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As given in Eq. (B-5) of Appendix B the gravity 
gradient torque for this orientation is 



gX \ 


T 

T 

- | 


gy 

l 


T i 

\ 


gz j 


T gmx s2 <* + *> c 
T gmy c { 3+4>) s2 ( n ) 

T gmz s (B+$)s2 


(5-1) 


where T , T , T_ represent the maximum values of torque 
gmx gmy gmz 

along the respective spacecraft principal axes and n is the angle 
between local vertical and the nodal axis , z 


n 


The torque 

components T gy and T gz are purely periodic whereas T gx is 
periodic about a bias value of* 


‘Vbias - V2 T gmx s2(6 + *> 


(5-2) 


Consequently, the magnitude of the gravity gradient torque- 
impulse (or bias momentum) over an orbit is* 

T 

I AH I = I / Tg dt l = l%^l|B2<B + i)| = AH max l s2 (B+$) I (5-3) 
0 


where T is the orbital period. This result applies for both 
circular and elliptical orbits. Over a mission |ah| will vary 
depending on 3. In Fig. (5-2) a typical variation of [ah|/AH = 

|s2(3+$) [ is shown for (j>=0° and 23° with the same orbit and 
launch conditions given previously in Fig. (4-1) . 

With momentum storage devices such as control moment 
gyros (CMGs) employed for attitude control, the periodic compo- 
nent of the gravity gradient torque-impulse is easily removed. 
However, the bias component of the torque- impulse causes an 
accumulation of angular momentum. In order to avoid eventual 
saturation, this bias momentum must be removed periodically. 
Three approaches for momentum dumping are based on: 


*This result is based on keeping 3 and ijj constant which is a 
reasonable approximation over an orbital period. 
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Pttfcg) 




AE WS VE SS AE 


AE — Autumn*! Equinox 
WS — Winter Solstice 
VE — Vernal Equinox 
SS - Summer Solstice 


BJItf+Sl $-23° 



Figure (5 - 2) - Typical Variation of j3 and Normalized Bias Momentum . 


1. reaction thrust 


9 

2 . magnetic control 


3. gravity field 10 ' 11 


large torques produced by 
reaction jets provide rapid 
desaturation ; 

interaction of a control magnetic 
dipole and the earth's magnetic 
field produces a bias torque for 
continuous gravity gradient bias 
torque cancellation or short 
period momentum desaturation; 

orbit dark- side maneuver and 
attitude hold operations produce 
gravity gradient torques which 
reverse angular momentum 
accumulation . 
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In the first two approaches | AH | given in Eq. (5-3) is 
the bias momentum dump requirement. In the third approach the 
dumping requirement will vary between |ah|/ 2 and | AH] , since 
the spacecraft is not held fixed in the solar-inertial mode 
during part of the orbit dark-side interval. For the first 
SVWS mission the gravity gradient dump procedure is superior. 
However, the other methods are useful alternatives when opera- 
ting, for example, with only two CMGs or during high 3 phases 
of a mission in which the dark-side time interval is small or 
non-existent.* Momentum dumping with magnetic control is also 
an attractive possibility for the second workshop and in 
future space stations where orbit dark-side maneuvers or reac- 
tion jet firings might upset critical telescope alignment 
during stellar astronomy experiments. 


5 . 2 Orbital Sunlight Interval 


The period during an orbit in which the sun is 
unocculted by the earth is referred to as the earth sunlight 
interval. In general, this interval depends on the proximity 
of the orbit path to the earth and its orientation with respect 
to the earth and sun line. Geometry for determining the shadow 
interval in circular orbits is shown in Fig. (5-3) where R is 
the earth's radius and H is the orbital altitude. Points P^ and 

T ?2 are locations in orbit where an observer's line-of -sight 

to the sun's center just becomes occulted by the earth. ** The 

earth shadow terminator corresponds to the great circle ABC which 
is normal to the sun line for a given 3. 

The orbit angle corresponding to the shadow interval 
is defined as 2a^. The corresponding angle for the sunlight ~ 

interval is 2 n where 

es 


"es = 180 ° - ®1 


(5-4) 


*Spacecraft maneuvers "during sunlight periods may interfere 
with solar experiments. As shown in Fig. (5-4) for circular 

orbits, the dark-side interval vanishes for 3>o=sin _1 (R/R+H) 
where R is the earth's radius and H the orbital altitude. 

**Umbra and penumbra effects of the earth-sun geometry are 
not included in this analysis. See footnote in Appendix E, 
p.E-3. 
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(b) 



Shadow 

Interval 


<2 <r,) 


Figure (5 - 3} • Geometry for Determining Earth Shadow Interval. 


as shown in Fig. (5- 3b) . The relation between and £ is 
obtained from a solution of the tetrahedron DEOP^ in Fig. (5-3a) . 

This yields* 


ccr^ = co/cB 


where 

so = R/(R+H) 

Hence, n_ is given by 
s s 

'cos ^(-ca/c£) 

n es ~ 

w 180° 


(B<cr) 


S< cr 


3>a 


(5-5) 


(5-6) 


(5-7) 


*In circular orbits co-^ is only a function of 6, which is 

assumed constant over an orbit. In elliptical orbits the shadow 
interval is a function of both 3 and i|i. Analytical results for 
this case are derived in Reference 13 and summarized in Appendix E. 
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Since the orbital angular velocity is constant in a 
circular orbit the fraction of an orbit in sunlight is 


T SL = ( 2 *es/ 360O) “ *es/ 180 ‘ 


(5-8) 


A plot of T gL as a function of 8 is shown in Fig. (5-4) for 

various orbital altitudes. When 8>a, the orbit is in total sun 
light. For a 230 NM circular orbit this occurs for 8>69.6°. 


50 8 ° 60 8 ° 66 9 ° 70 9 ° 76 3 ° 


Figure {5 - 4) - Fraction of a Circular Orbit in Sunlight vs^/3 
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5.3 Electrical Energy from Solar Arrays 

Electrical energy from solar arrays is a function of 
the array orientation relative to the sun line and the duration 
of the sunlight interval. In a solar orientation, such as the 
solar- inertial mode illustrated in Fig. (1-1) , the incident solar 
energy on the arrays is a maximum. In orbit-oriented modes, 
such as LV, LH and POP,* the maximum is obtainable only with an 
independent array pointing capability. This may be achieved by 
two-degree-of- freedom (TDOF) articulation of the array or by 
combining single-degree-of- freedom (SDOF) array articulation 
with spacecraft rotation about the roll axis. However, contin- 
uous articulation of large solar arrays at orbital rates in 
long duration missions may be undesirable from a reliability 
standpoint. Nevertheless, an improvement over fixed arrays can 
still be achieved with SDOF arrays which need only be articulated 
intermittently to account for the relatively slow (daily) varia- 
tion in g. 

In the following subsections the effect of p on the 
electrical energy available from fixed and incrementally articu- 
lated (SDOF) solar arrays is examined for LV, LH and POP attitude 
modes and various spacecraft roll profiles. In each case optimum 
panel articulation for maximum utilization of the incident solar 
energy is evaluated as a function of 0. The subsequent analysis 
is based on the following assumptions: 

(1) secondary effects on array output such as fluctua- 
tions in panel temperature and loss of incident 

14 

energy due to cover glass reflection are ignored, 

(2) the solar panels are mounted such that shadowing 
by the spacecraft is negligible, 

(3) the panel articulation axis is normal to the space- 
craft roll axis , and 

(4) the solar panels may be rotated only intermittently 
(e.g. daily). 


*L V, LH and POP modes are defined in Footnote **, p. 1. 
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5.3.1 Evaluation of Electrical Energy 


The electrical power (P) generated by a planar solar 
array for an arbitrary array orientation is expressed as a 
fraction of the maximum electrical power (P^) which is generated 

when the array is pointed directly at the sun. This relationship 
is given by 

P = P M cX |X|<_90° (5-9) 

where cA is the cosine of the angle between the sun line and the 
outward normal to the active side of the solar array. Generally, 

cA is a function of four angles: , r»/ <{> t and a . As illustrated 

P • 

in Fig. (5-5) , n is the current location in orbit relative to noon, 

<f> is the spacecraft roll angle and a is the array articulation 

P 

angle. * 



Figure (5 - 5) - Spacecraft Location and Solar Array Orientation Angles. 


*The reference for a =0° 

P 

ly for each attitude mode. 


and <J>=0° will be defined subsequent- 
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The array power output is zero whenever the space- 
craft is within the earth shadow interval ln|>n * or whenever 

1 1 es 

cX<0, which represents panel self -shadowing. The respective 
points in the orbit ( r > ==ir i SS 2 r r,=n ssl^ corresponding to entry into 

and exit from the self -shadow region can be obtained, by solving 
the equation, cX=0, for n . Array self-shadowing is of interest 
whenever r i ss 2 <n es a nd/or n ss ^ >_ n es . ** In that event the self— 

shadow region includes part of the sunlight region, as shown in 
Fig. (5-6) for two possible cases, and thus reduces the total 
electrical energy obtainable from the solar array. 


*1 - 





V = V SS 2 



l b > ss2^ | *?$sl| 


Legend 

m Earth Shadow Region 
tSS>j Panel Self Shadow Region 


Figure (5 - 6) - Earth Shadow and Array Self-Shadow Regions in Circular Orbits 


The electrical energy (E) available per orbit from 
the solar array is equal to the integral of the output power 
(P) over the period in the orbit when the panels are exposed to 
sunlight. This integral can be written as 


E 



(5-10) 


where correspond to the point of array entry into shadow, 

<n r ti) corresponds to the point of exit from shadow and n is 


*See Section 5-2 and Fig. (5-3) . 

**In general, ln ss2 l ^ |n ssl |/ since the solutions of cX=0 

for n are not necessarily symmetric about n=0°. However, for 
many cases considered in subsequent work, n ss 2 = “ n ssl Hn ss* 
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the instantaneous angular velocity in orbit. The average power 
supplied during the period C^-t^) is simply 


= 1 2pat = 4*1 (5 - U) 

t l 

For circular orbits, is constant and = ( T) 2 ~ n l) ' 

so that the expressions for the electrical energy and average 
power become 



and 

n 2 

cAdn 

1 


OJ 


ave 


E = 


-M 


(n 2 


I 


(n 2 -T ii) "n 


(5-12) 


(5-13) 


where Eq. (5-9) has been substituted into Eq. (5-12) for P. 

The maximum electrical energy (E ) is obtained with 

s 

the panels pointed continuously at the sun, which means that 
cA=l. Since panel self-shadowing is precluded in that case, 
ru=ri _ , t)-, = -n and Eq. (5-12) yields 

Cm GS X 05 


p p 

w _ M , . M , 0 . 

E s = — ^2~^1 = — 2tl es 

o o 


(5-14) 


Hence, the maximum electrical energy obtainable from solar 
arrays varies with $ and orbital altitude as shown by the curves 
in Fig. (5-4) for n . 

GS 


For an arbitrary panel orientation the available 
electrical energy given by Eq. (5-12) can be expressed in terms 
of the sun-oriented array energy output (Eg) by substituting for 

P M /t» 0 from Eq.'(5-14) . This yields 




cAdn 


E 


(5-15) 
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or in normalized form 


E = (E/E s ) * 


2n 


es 


r 


cXdn 


(5-16) 


Thus , E represents the fraction of E g obtainable for a given 
attitude mode and array orientation.* 


In view of Eq. (5-13) the maximum average power 
obtainable with the array pointed continuously at the sun is 
just P . Consequently the normalized average power is given by 


r n 2 

p = (P /p ) = — - — r- cXdn = 

v ave / m (ru-ni ) 


2n 

es 


(5-17) 


This indicates that the normalized electrical energy , E and the 

A 

normalized average power P are identical, if self— shadowing 
occurs only during the earth shadow region (i.e., r, 2"" n l =s ^ r, es^ * 


*An alternate normalization of the array energy output (E) 
could be made relative to the sun-oriented array output (E go ) at 

3=0°. In view of Eq. (5-14) it follows that E so = ( p m / w 0 ) 2n eso 

where 2n corresponds to the sunlight interval at 3=0° • 

©SO 

Accordingly n in Eqs. (5-15) and (5-16) is replaced by n and 
the normalized array output is designated E q . Also, if E g in 
Eq. (5-14) is normalized to E gQ , the sun-oriented array output 
can be expressed as 


E he /E 
s s' 


so 


=n 


es 


/n 


eso 


which varies with 3 analogous to the curves for n es /l80° in 

Fig. (4-4) . Normalization relative to E g is preferred in this 

section however, since this approach will indicate directly the 
array output relative to the sun— oriented array for all 3» 

Nevertheless, the two approaches are simply related, since 

A A A 

E =E E. 

O S 
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5.3.2 Electrical Energy in Local Vertical Modes 

The solar array orientation in the LV mode is illustrated 
in Fig. (5-7), where the panel coordinate axes, y and z , define 

Jr Jr 

the articulation axis and the normal to the active side of the 
array. When the panel angle (a^) and the spacecraft roll angle 

(<j>) are zero, the array lies in the orbital plane with x parallel 

ir 

to the spacecraft roll axis (Xp Q ) . For arbitrary and (j> the 
expression for cA given in Eq. (D-19) of Appendix D is 

cA^y = Cct p c< f* + sa p Ct 1 ~ Ca p s< f > sri (5—18) 

In the following subsections E is evaluated as a function of 8 
from Eqs. (5-16) and (5-18) for various panel angles and several 
spacecraft roll profiles.** 



Figure (5 - 7) - Solar Array Orientation in Local Vertical Mode . . 


*The results in this section were contributed previously to 
Reference 15. 


**'As may be surmised from Fig. (5-7) it is really only 
necessary to consider 0°<6<90°, since redefining the spacecraft 
roll such that <j>-*-180°-<|> , when -90°<8<0°, yields exactly the same 
expression for cA LV in Eq. (5-18) . Also, in view of Fig. (5-7) it 


is sufficient to consider 0 o £ctp<90°, since in the range, 

-90 0 <_otp<0° , serves only to orient the outward normal’ to the array 

(z ) further away from the sun line. 

P 
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5 . 3 . 2 . 1 Solar Arrays Operated with Spacecraft Roll (I/O ) 

The angles and $ corresponding to a solar oriented 
array (cA^-l) are defined by Eqs. (D-ll) and (D-12) , namely 

sa = cn c$ (5-19) 

XT 

tamj> = -sri/tan$ (5-20) 

Although variation of at orbital rates is precluded by 

assumption, the variation of <j> according to Eq. (5-20) represents 
the optimum spacecraft roll in the sense that this maximizes the 
instantaneous array power output.* A plot of $ ^ as a function 

of n is illustrated by the solid curves in Pig. (5- 8a) for various 
values of 6 . 


(a) (b) 




Figure (5 - 8) - Spacecraft Roll Angle vs Orbital Position Angle for Three Roll Modes. 


Although theoretically possible, the variation of 
^ according to Eq. (5-20) imposes some adverse conditions on the 
spacecraft attitude control system. One practical aspect is that 
the spacecraft must follow a rather complex roll command, which 
varies with $. The primary disadvantage, however, is that the 
spacecraft must execute an abrupt roll maneuver twice per orbit, 
when g is small . * 


*Since P=P m cA, it follows that 3/3<}> (cX LV ) =0 yields Eq. (5-20). 

*In view of Figure (5- 8a) the theoretical maximum roll rate, 

I I max ' occurs at n=0° and ±180° approaches infinity as 

8-*-0°, since ||| =to /tan 8. 

1 T 'max o 
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Two easily implemented suboptimal alternatives to the 
optimum roll are represented by the dashed curves in Fig. (5-8a) . 
In one case the roll motion (<J> cont ) is continuous with 

. In the other case the roll motion (<f> ) is cyclic 

conu o cy c " A 

about zero such that 


n+180° -180° < n < - 90° 

-n - 90° ± n <_ 90° 

n-180° 90° < n < 180° 


(5-21) 


As is evident from Fig. (5-8a) both cases closely approximate 
<J>opt near n=0° and 8=45°. 

These two alternatives will be compared with the 

optimum roll case in terms of the electrical energy (E) available 
from the array as expressed in Eq. (5-15) . The boundaries of the 

a 

array sunlight interval needed for evaluating E, can be 

determined by comparing the earth shadow boundary ( ± n _ ) with the 

W to 

panel self-shadow boundary (n ssl rn ss 2 ) described previously. In 

Appendix C the self-shadow boundary is shown to be symmetrical 
about n-0° for all three cases so that n -n =n . Conse- 

SS<d SSI ss 

quently the array sunlight interval is determined by 



The relationship between n _ and n _ for various values of a 

os s s p 

and 8 is shown by the curves in Fig. (C-l) , (C-2) and (C-3) in 


•Appendix C for the optimum roll, continuous roll and cyclic roll 
cases respectively. 


With these data the normalized electrical energy (E) 
can be evaluated after integration of Eq. (5-16) with Eq. (.5-18) 

A 

substituted for cA. The normalized average power (P) available 
during the sunlight interval is then obtained by using Eq. (5-17) . 

A 

The expressions for E as given in Eqs.(C-12), (C-19) and (C-30) 

are listed below. 



BELLCOMM, INC. 


48 


Optimum Roll * 


1 f 

= < sa 

n es l 

= (' 1 es /n SJ ) E 


p c$ sn SJl + cotp [E (k) + E (k. 


0)]j 


where 


P = 


k = cB 


and 


9 - n s C 90 ° 
Continuous Roll 


E = 


1 r ca eg 

U ) — 2s (a p +e)sT1 s)l + — ’ — <2n„.-s2n..)| 

es L 


‘SA SA 


i " (W n s*)= 

Cyclic Roll 


E = 




ca cB 

4ca p sB + 2s(a p >B)sn s£ + — (2n s £ "s2n SJl ) 


(5-23) 

(5-24) 

(5-25) 

(5-26) 

(5-27) 

¥ 

(5-28) 

(5-29) 



(5-30) 


In Pigs. (5-9) , (5-10) and (5—11 ) , E is plotted as a 

function of B and the panel angle (a ) for the various roll 

ir 

modes. The corresponding results for P differ significantly 
from those of E only in the continuous roll mode. The curves of 
P vs B for this case are shown in Fig. (5-12) . The knee in the 

A 

curves of E vs 3 near 3=70° is due to the rapid decrease in the 
earth shadow interval.** In each case the upper envelope boundary 

A 

on the curves for E represents the maximum electrical energy 

A 

(E max ) available from the array with optimum variation of with 

3. A comparison of optimum results for these cases and optimum 
results in other attitude modes are given in Section 5.3.5. 


*The function E(k/0) in Eq. (5-23) is an elliptic integral of 

16 

the second kind which is available in tabulated form as a func- 
tion of the modulus k and argument 0 . The function E (k) is the 
complete elliptic integral of the second kind (0=90°) . 


**See Fig. (5-4) . 
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/3{deg) 

Figure (5 - 1 1) - Normalized Electrical Energy From Solar Array as a Function of a p and /3 
(LV - Continuous Roll Mode). 
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0{Deg) ' 

Figure (5 - 12) - Normalized Average Power from Solar Array as a 

Function of a p and j3 (LV - Continuous Roll Mode) 
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5 . 3 . 2 . 2 Solar Arrays Operated with. Fixed Spacecraft Roll ($=0 ) 

If the spacecraft is roll stabilized, the roll degree- 
of-freedom for pointing the solar array is not available on a 
continuous basis. Unless otherwise constrained however, the 
spacecraft roll angle is arbitrary and can be fixed at some 
desired value ( <f>=constant) . * It is of interest therefore to 

a 

evaluate the electrical energy (E) available from the array for 
various fixed values of <(» and to determine as a function of $ 
the optimum combination of fixed roll angle (<(>) and panel angle 

A 

(a ) which maximizes E. 


As in previous work the sunlight interval must be 

A 

determined in order to evaluate E according to Eq. (5-16) . In the 
fixed roll mode the panel seif-shadow boundaries (n s n S s2^ are 

not symmetrical about n=0° except when $=0 9 . Consequently, the 
self-shadow and earth shadow intervals may assume any one of four 

possible arrangements, as shown in Fig. (5-13) . 


A) 


C) 


B) 


17 = 0 (Noon) 



► Sun 



Orbit Path 



► Sun 


ss 2 ^ ssl Orbit Path 


Legend 


I Earth Shadow Region 
(S3 Panel Self-Shadow Region 


^ Sun 

V= 0 
(Noon) 


Figure (5 - 13) - Possible Earth Shadow and Solar Panel Self-Shadow Regions (LV - Fixed Roll Mode). 


*If the spacecraft is roll stabilized by reaction thrust 
methods, <j> will normally drift within a limit cycle deadband. If 
the deadband' is sufficiently small however, <}> can be considered 
essentially constant. If momentum exchange devices are used for 
control, a precisely fixed value of <j> can be maintained. 
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In general the self-shadow boundaries (n .,11 o) vary 

with 3. In Figs.(C-5), (C-6a) and CC-6b) of Appendix C.4 the vari- 

ation of n ss -^ and n ss2 as a function of 3 and is shown for 

<j>=-45 o ,0° and -90° respectively. Also shown are the earth shadow 
boundaries ( * n ) corresponding to a 230 NM circular orbit. The 

various shaded regions (A, B, C or D) in Figs.(C-5) and (C-6) 
correspond to the respective shadow boundary arrangements in 
Fig. (5-13). 

A 

With this information the electrical energy (E) can be 
evaluated after integration of Eq. (5-16) . In order to include 
the effect of distinct earth shadow and panel self-shadow inter- 
vals as in Fig.(5-13d), it is necessary to modify this expression 
so that 


E 


1 

r n 2 

| I c;>l LV^ n 

r n 2 

CX LV dn j 

► (5-31) 


2 n ^ 

es 

" J. 





J 


1 I 

the limits (n^n-^; n 2 / r i 1 ) 

corresponding 

to the various 

in Fig. 

(5- 

•13) are defined 

as 





f n ss2 


Case: 

A 

(5-32a) 

n 2 


n es 

V 


Case: 

B , C,D 




r ^ssl 


Case: 

A, C 

(5-32b) 

n l 


J" n es 


Case : 

B,D 


1 


r^ssi 


Case: 

D 

(5-33a) 

n 2 


1 

[ 0° 


Case: 

A,B,C 


1 


r r| ss2" 

-360° 

Case: 

D 

(5-33b) 

n l 

H 

l 

0° 

V. 


Case: 

A,B,C 
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Substitution for cA from Eq. (5-18) and integration of Eq. (5-31) 
yields 

E = — [A(n 2 ~n 1 ) + BCsn^-sr^) + C(cri 2 -cn 1 )] 

n es 

- [A ( n 2 _ ^i^ + 6 (sn 2 — sn 2 ) + 0 ( 0112 - 0112 ] (5-34) 

where 


A = cot c<j> sB 
P 

(5-35) 

B = sa cB 
P 

(5-36 ) 

C = +ca cB s<j> 
P 

(5-37) 


In Fig. (5-15a,b, c) E is plotted as a function of 6 and the panel 
angle (a p ) for three different roll angles: <j>=0°*, <{•= -45° and 

4 .= -90°. 


Three particular cases are of interest, since they 
represent possible orientations of a spacecraft in a local verti- 
cal storage mode. The three panel orientations and the associated 
<Xp and <J> are shown in Fig. (5-14) . The corresponding shadow bound- 
aries for each case are given in Fig. (C- 6 ) in Appendix C.4. 

A 

Consequently, the respective expressions for E obtained from 
Eq. (5-34) are as follows: 


1 ) g p =o°, H° (n 2 = ~ r| 1 ~T^ es ) ! 


E = 2n ~ [s g ( n 2 ~ n x ) ] = s3 


(5-38) 


2) ot= 0 0 , <j>= -90° (n 2 = + n es /n 1 =0°): 


E = 


2n 


es 


[- 03 ( 0112 - 0112 )] = (c3/n es )(l-cn )/2 (5-39) 


*This case was studied in Reference 17. 
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3) a p =90 °/ H° Cn. 2 - -11^90°) : 


E 



[eg (s tt_ 2 s ) 1 


= cg/n 


es 


(5-40) 


a) a p =■ 0°, 4> = 0° 


c) a p - 90° <t> - 0° 


Local Vertical Loca1 Vertical 



Figure (5 - 14) - Three Possible Panel Orientations for Spacecraft in a Local Vertical Storage Mode. 

The curves in Fig. (5-15) indicate the electrical 

A 

energy (E) available for other panel orientations as a function 
of g. The upper envelope boundary on each family of curves 

A 

represents the maximum electrical energy (E ) available from 
the array with optimum variation of with g. Plots of E raax an< ^ 

vs g obtained numerically for several fixed roll angles (<f>) 
are shown in Figs. (5-16) and (5-17) respectively. These results 

* A 

indicate that <j>=0° apparently yields the largest E max for all g. 

As discussed in Appendix C.4, however, other combinations of a 

^ IT 

and <j> exist, which yield the same E max as with <j>=0° for g in the 

range 0°j<g<23 o . Bounds on a ^ for other possible values of $ 

are shown in Fig. (C-8) in Appendix C.4. From an attitude control 
standpoint however, the simplest scheme would be to use <|> t =0° 

for all g . 
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Figure (5 - 16) - Maximum Electrical Energy from Solar Array as a 

Function of ( 3 . (a p = a p ; LV - Fixed Roll Mode). 



Figure (5 - 17) - Optimum Panel Angle vs- jS (LV - Fixed Roll Mode) . 
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5.3.3 Electrical Energy in Ixpcal Horizontal Modes 

The solar array orientation in the LH mode is illustrated 

in Fig. (5-18) with the same panel coordinates (x ,y ,z ) as used 

P P 3? 

in the LV mode. When the panel angle (a^) and the spacecraft roll 

angle ($) are zero, the array lies in the orbital plane with y , 

P 

the articulation axis, directed along the upward local vertical 
(z^) . For arbitrary and <f> the expression for cX given in 

Eq. (D-20) of Appendix D is 



Figure (5 - 18) - Solar Array Orientation in Local Horizontal Mode 


In the following subsections E is evaluated as a function of 6 
from Eqs.(5-16) and (5-41) for various panel angles and several 
spacecraft roll profiles.* 


*As in the LV mode it is only necessary to consider 


0<6<90° and 


0<a <90° 

- p- 


See Footnote**, p. 45 . 
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5 • 3 . 3 . 1 Solar Arrays Operated With. Spacecraft Roll ($/0 ) 

The angles and <{> corresponding to a solar oriented 
array (cA LH =l) are defined by Eqs. CD-13) and (D-14) , namely, 

sa = -sn eg (.5-42) 

P 

tan <j> - -cn/tan g (5-43) 

Although is restricted to incremental variation, spacecraft 

roll according to Eq. (5-43) is optimum in the sense of maximizing 
instantaneous array power output.* A plot of as a function 

of n is illustrated by the solid curves in Fig. (5-19) for various g. 



Figure (5 - 19) - Spacecraft Roll Angle vs- Orbital Position Angle for Two Roll Modes. 


*Analogous to the local vertical case 3/3 <j> (cX LH ) =0 yields 
Eq. (5-43) . 
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The variation of <p according Eq. (5-43) leads to abrupt 
roll maneuvers twice per orbit when g is small, just as in the 
LV mode. A practical suboptimal alternative* to optimum roll is 
the cyclic roll profile shown in Pig. (5-19) and defined by 


4 > 


eye 


'-90°-n 

-90°+n 


-180°<n<0° 

0 o <n<180 o 


(5-44) 


This approximates ^ closely near n= ± 90 l 


and 8=45°. 


As found in Appendices C.5 and C.6 for the LH optimum 
and cyclic roll modes, the panel self-shadow boundaries 
^ T1 ssl' T1 ss2^ are s Y jnmetr ^ cai about n=+90° rather than n=180°. 

Consequently Eq. (5-16) for evaluating E must be modified slightly 
to accommodate the possibility of distinct earth shadow and self- 
shadow regions as illustrated in Pig. (5-20) . Thus, 


E = 


2ti 


es 


r n es r " 

J CX LH dn “ J 


c\ LH d n 


-ri 

v es 


where 




n 2 " 


es 


ss2 


^es— n ssl 


l T1 s S r (180 ' > ~ n ss2 ) n es >1 ? S sl 


(5-45) 


The relationship between n and n o 

BS 5S» 


(Case C) 

(Case D) 

for various values of a 


(5-46) 


and 8 is shown by the curves in Figs. (C-9) and (C-10) for the 
optimum roll and cyclic roll cases. 


*In the LH mode a continuous roll profile analogous to that 
used for the LV mode does not yield a sufficiently close approxi- 
mation to * . and is therefore not considered, 

opt 
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*? = 90° 


V= 90° 




Legend 

/////' Earth Shadow Region 

Panel Self Shadow Region 


(A) ^es < ^ssl (180° - 


(8) V > ’Jssl = ( 180 ° - »fe2> 


Figure (5 - 20) - Possible Earth Shadow and Panel Self-Shadow Regions in 
Local Horizontal Optimum and Cyclic Roll 'Modes. 


The corresponding evaluation of E for cases C and D as given by 
Eqs.(C-93) and (C-100) of Appendix C yields 

Optimum Roll * 


Cl E = [ca p [E(k,n es ) + E(k,n ss2 )l + s« p ceCcn ss2 -c,, es )j 


D) E 




(5-4 7) 


Od p IE(k) + E(k,„ es ) - E(k,n ss2 ) ] + s« p c6cn ss2 


where k=cB , 

Cyclic Roll 


C) E 


= 1 f 

2 ^esV 


t 2ca p s 6 - s(a p+ B)cn es + s(a p -6)cn ss2 


) 


Co cs 


[2 (r 'es +r) ss2 ) + (s2r ’i=B +s2r i <3a 9) ] 


es n ss2' 


■) 


D) E 


“ + st V e)on ss2 + 


ca eg 

+ -^—[2 (n es +n ss2 -ir/2) + (s2n es +s2n gs2 ) ]J 


(5-48) 


*See Footnote *, p.48 
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>3 (deg) 

Figure (5-21) - Normalized Electrical Energy from Solar Array as a Function of Op and j3 
(LH - Optimum Roll Mode), 


Sun Oriented Array E - 1 0 



(deg) 


Figure (5 - 22) - Normalized Electrical Energy from Solar Array as a Function of a and j3 
(LH - Cyclic Roll Mode) P 
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In Figs. (5-21) and (5-22), E is plotted as a function of 3 and 
a for the optimum roll and cyclic roll modes. In each case 

the upper envelope boundary on the curves for E represents the 

A 

maximum electrical energy (E ) available from the array with 

max 

optimum variation of a with 3. A comparison of the optimum 

ir 

results is given in Section 5.3.5. 

5 . 3 . 3 . 2 Solar Arrays Operated with Fixed Spacecraft Roll (j>=0 ) 

The local horizontal fixed roll mode is of importance 
because it is suitable for earth pointing experiments , it is 
easily acquired from the solar-inertial mode and external 
disturbance torques acting on the spacecraft tend to be small. 
It is of interest therefore to evaluate the electrical energy 

A 

(E) available from a solar array for various panel and fixed 
roll angles (otp,(f>) and to determine a combination which maxi- 
mizes E. 


As in previous work the array sunlight interval must 

A 

be determined in order to evaluate E. Analogous to the local 
vertical fixed roll mode the array self-shadow boundaries 
(n ss i, r i SS 2 ^ are not symmetrical about n=180°. Consequently, 

the self-shadow and earth shadow interval may assume any one 
of four possible arrangements shown in Fig. (5-23) . In 
Figs. (C-ll) , (C-12a) and (C-12b) of Appendix C the variation 


•45* 


0° and -90° 


of n , and n 0 with 3 and a is shown for <j>= 
ssi ssz p 

respectively.* Also shown' are the earth shadow boundaries 

corresponding to a 230 NM orbit. The various shaded regions 

(A,B,C or D) in Figs. (C-ll) and (C-12) correspond to the 

respective shadow boundary arrangements in Fig. (5-23) . 


*For a given and 6 the self-shadow boundaries in the 

LH mode are rotated 90° back (clockwise) toward noon (ri=0°) 
from corresponding boundaries in the LV mode. This follows 
because substituting Ti=n-90° into Eq. (5-41) for cX LH yields 

c\ LV in Eq. (5-18) . 
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Legend 

\///X Earth Shadow Region 

Panel Self-Shadow Region 
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ssl 

Case: 

D 

(5— 50c) 

0° 

Case : 

A,B,C 


ss2 

Case: 

D 

(5-50d) 

0° 

Case: 

A,B,C 



Substitutions for cX LH from Eq. (5-41) and integration of 
Eg. (5-49) yields 


A 

E 



IA(n 2 -n 1 ) + Btcr^-cn^) 


II t I 

”* (A ^ ( cti 2™ dl ) — C 


- C (sn 2 _sr i 1 ) J 

1 1 

(sr^-sn-^) ] 


(5-51) 


where A, B and C are as defined in Eqs. (5-35) , (5-36) and 
(5-37) for the local vertical mode. In Fig. (5-24a,b,c) E is 
plotted vs 8 for various panel angles (a p ) and three different 

roll angles: <f>= 0°, -45° and <f>= -90°. The upper envelope 

boundary on each family of curves represents the maximum 
electrical energy (E raax ) available from the array with optimum 

variation of a with 8. Plots of E and a , .. vs 8 for 

p max p (opt) 

several fixed roll angles are shown in Figs. (5-25) and (5-26). 

The upper boundary on the curves in Fig. (5-25) 
represents the largest E obtainable with optimum variation 

IuQA 

of both cip and <j>. As discussed in Appendix C.7 this can be 

achieved with a p =0° and $=* -<* ( for the LV-Fixed Roll mode 

as shown in Fig. (5-26) . As in the local vertical case however, 
the optimum combination of a and <{> is not unique for O<0<23°. 

P A. ^ 

Bounds on a , .. and <j> . which yield the same E are shown 

p v.opt; opt max 

in Fig. (C-13) . From an array design standpoint the simplest 

scheme would be to use a fixed array (a =0°). 

P 
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Figure (5 - 25) - Maximum Electrical Energy from Solar Array as a 

Function of 0 < a p = a p { 0 p t )' LH- Flxed Ro11 Mode ) • 



Figure (5 - 26) - Optimum Panel Angle vs /3 (LH - Fixed Roll Mode)- 
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5.3.4 Electrical Energy in Perpendicular-to-Orbital-Plane Modes 

The solar array orientation in the POP mode is illu- 
strated in Fig. (5-27) with the panel coordinates (Xp,y , z ) 

defined as in the LV and POP modes. When the panel angle (a ) 

P 

and the spacecraft roll angle (<j>) are zero, the array is normal 
to the orbital plane with directed’' along the upward local 



Figure (5 - 27) - Solar Array Orientation ih POP Mode 


For arbitrary a p and 4> the expression for cA given 
in Eq. (D-21) of Appendix D is 

cXpOp = c 0 c(<|>+n) - sctp sB (5-52) 

Two POP modes of interest are based on a solar oriented array 
(S/POP) and a local vertical oriented spacecraft -(LV/POP) for 
accommodating earth pointing experiments. 

S/POP Mode 

The angles a and <}> corresponding to a solar oriented 
P 

array (cAp 0 p=l) are simply 

a = - B (5-53) 

p 


<j> = -n 


(5-54) 
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so that E in Eq. (5-16) is just E=1.0. This mode was to be used 
for AAP-1/AAP-2/AAP-3A missions with the "Wet Workshop" config- 
uration . 

LV/POP Mode 

In a LV/POP mode <f>=0° is the spacecraft roll angle 
for accommodating earth pointing experiments with a line-of- 
sight normal to the array articulation axis (y p ) . Hence, 

cA pnp = ca eg cn - so sg = so eg cri + ca sg (5-55) 

jtkjst p P ir P 

which is equivalent to cA LV in Eq. (5-18) for <j>=0° r if the 

angle a =90°+a shown in Fig. (5-27) is interpreted as the 

P P _ 

panel angle. Thus, all results in Section 5. 3.2.2 for the LV- 
Fixed Roll mode (<f»=0°) are applicable to the LV/POP mode. 
Specifically this includes curves for E, E max and vs 8 

given in Figs . (5-15a) , (5-16) and (5-17). 


5.3.5 Summary and Comparison of Optimum Solar Array Performance 

Solar array energy output (E) normalized to the output 
from a sun-oriented array at the same g has been evaluated for 
LV, LH and POP spacecraft attitudes and various roll profiles. 

In each case E is formulated as A a function of g and the panel 
angle (a ) . For 230 NM orbits E is plotted vs g for various 

values of a in Figs. (5-9,10,11,15,21,22 and 24). The upper 
P 

envelope boundary on these curves represents the maximum electrical 
energy (E ) available from the array with optimum variation of 

Itl&X A 

a with g. Plots of E and the corresponding optimum panel 
p max 

angle and roll angle (fixed roll modes) are shown in Figs. (5-28)- 
(5-31). Analytical expressions for E max , a p (opt) and (}> opt (fixed 

roll modes) are given in Appendix C. For comparison, the maximum 

electrical energy (E ) normalized to the sun-oriented array 

o max 

output at 8=0° * is shown in Fig. (5-32) . 


Solar array performance in the LV-Optimum Roll mode is 
up to 50% higher than in the LH-Optimum Roll mode and at worst 
is only 15% less than the sun-oriented array performance. In 
the LV and LH-Cyclic Roll modes performance approaches 85% of 
that for the sun-oriented array in the intermediate range of 8 
(30°<8<60° ) with the LV-Cyclic Roll mode significantly better than 
the LH-Cyclic Roll mode for 8 <30°. 


*See Footnote *, p. 44. 
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In fixed roll cases optimum array performance is the 
same in all three attitude modes: LV, LH and POP. The results 

indicate that the maximum array output is only 52% of the sun- 
oriented array output at 0=0°, but increases steadily with 
increasing 0. Maximum output is achieved in LV and POP modes 
with array articulation only U O p t =0°/ a p (opt) = ^opt^ * w ^ ereas 

in the LH mode, only incremental spacecraft roll is necessary 

( *opt = -So P t' Voptr° 0> - 

These results represent the maximum which can be 

achieved with array pointing using only spacecraft roll and/or 

incremental SDOF array articulation,** The curves of E vs 0 

u max 

permit a direct assessment of the minimum possible performance 

degradation in several orbit-oriented modes from the performance 

of solar-oriented arrays.*** Additional degradation for 

particular configurations where spacecraft shadowing of arrays 

is a significant factor must be evaluated separately. Changes 

in orbital altitude of perhaps ±30 NM from the 230 NM altitude 

used in this analysis will tend to have only slight effect. 

Altitude affects n , which in view of Fig. (5-4) varies only 

©S " ^ 

slightly with a ±30 NM altitude change. For larger altitude 

changes the analytical results in previous sections and in 

Appendix C can be used to evaluate E and Ei vs 0. 

max 


*See Appendices C.4 and C.7 for derivation of £ 0 p t and 
evaluation of * opt and « p(optK 

**Results in Section 4 for |p| max and I £ l max ma Y be used 
with curves of a ^ vs 0 in assessing maximum articulation and 

articulation rate requirements. 


***The solar array model used in this analysis, Eq. (5-9) , in- 
cludes the primary effect (cX) on array output. When the sun 
incidence angle (X) is large (e.g. X>70°) , temperature fluctuations 

and cover glass reflection^ may have a significant effect on the 
instantaneous power output. However, unless X is large over an 
appreciable portion of the array sunlight interval, the impact of 
these effects on energy output per orbit is relatively small. 
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Figure (5 - 28) - Maximum Electrical Energy from Solar Array vs (3 
(LV and POP Modes). 



Figure (5 - 29) - Maximum Electrical Energy from Solar Array vs /J 
(LH Modes). 











o max' 
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Figure (5 - 32) - Maximum Electrical Energy from Solar 
Array vs. J3 (Normalized to Sun-Oriented 
Array Output at j3 = 0°). 


1.63 
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5.4 Incident Thermal Energy on Spacecraft 


Thermal energy - incident on spacecraft surfaces in 
earth orbit is derived primarily from three external sources: 
direct solar radiation, earth reflected solar radiation and 
earth-emitted thermal radiation. In solar-oriented attitudes 
the direct solar component is constant during sunlight intervals 
whereas the reflected solar and earth thermal components vary 
with location in orbit and B. In earth-oriented attitudes 
the direct and reflected solar components vary with location 
in orbit and B while the earth thermal component is constant 
for circular orbits and varies only with altitude for ellipti- 
cal orbits . 

The incident radiation from each energy source can 
be evaluated in terms of view factors which account for the 
surface configuration and other parameters such as orbit 
altitude and the orientation of the surface with respect to 
the sun and the earth. View factors are normally derived on 
the basis of the following assumptions : 

1. direct solar radiation impinges with parallel 
rays, 

2. earth emitted radiation is diffuse and con- 
stant for all points on the earth (earth 
assumed in state of thermal equilibrium) , 

'3. earth is a diffuse reflector with constant 
albedo over the entire surface, and 

4. atmospheric scattering effects of direct 
solar radiation on spacecraft surfaces 
are negligible. 

In terms of view factors the incident thermal radia- 
tion from each energy source can be expressed as 


Q d SA‘F ds « 

(Direct Solar) - 

(5-56) 

Qr = sa#1f rs 

(Reflected Solar) 

(5-57) 

Grt sa - T~ f et 

(Earth Thermal) 

(5-58) 


2 12 

Here S represents the solar constant (443 Btu/hr/ft ) , A is 

the maximum projected surface and a is the earth albedo con- 
12 

stant (0.36)’. The various view factors are represented by 



Bellcomm, Inc. 


74 - 


F_ 0 , F n „ and F_ m .* Since the spectral content of earth thermal 

Ub Kb i li 1 

radiation is wide, it normally is treated separately whereas 
direct and reflected solar radiation frequently are combined. ** 

The incident thermal energy per orbit upon a surface 
is obtained by integrating each component of the incident 
thermal radiation over an orbital period (T) : 


W. = 
1 


i*T/2 

I Qj_dt - I Q i dn/n 

-T/2 Iff 


(i = DS ,RS or ET) 


(5-59) 

where n is the instantaneous orbital angular velocity. In the 
following work only circular orbits are considered so that 
n = o o =2ir/T. It is convenient to normalize with respect to 

the direct solar component for 8=0° with the area A normal to 


the sun line (i.e., 

— n < n < n ) . *** 

eso eso 


F ds = 1.0 during the sunlight interval 
In normalized form 



V W n 



(i = DS,RS or ET) 

(5-60) 


*As a consequence of the parallel ray assumption for 
direct solar radiation, F DS is merely the direction cosine 


between the sun line and the positive normal to area A. Evalu- 
ation of F^g and F ET entails integration over the surface of 

the earth which can reflect or radiate thermal energy to area 

A. Heretofore, F G and F™ have been evaluated numeric- 

18 19 20 18 

ally ' 1 and tabulated for certain ranges of parameters. 

21 

Recent analytical results have led to closed form expressions 
which cover several geometrical shapes and a wide range of 
parameters exactly and approximate well the exact case in other 
ranges. 


**An equivalent view factor corresponding to the sum of 
direct and reflected solar radiation is given by 

F S = F DS + ^RS s V SA 


***Here ±n represents the orbital terminator at 6=0°. 

eso , 

See Fig. (5-3*b) . 
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where 

W n * TP f QdS 31 ' ” JT 
o y c 


Specification of a particular attitude mode and 
surface configuration leads to particular expressions for the 
view factors F Dg , F RS and F ET which are functions of 8 , orbital 

position ti and orbital altitude H.^ Integration of Eg. (5-60) 
for each radiation component yields the normalized incident 
thermal energy. 

5.4.1 Cylindrical Surface Configurations in Earth-Oriented Modes 


( SS ° Sad " - !f( 2 "eso) < 5 - 61 > 

_n eso 


The incident thermal energy is evaluated first for 
LV, LH and POP attitude modes* and a cylindrical spacecraft 
model with hemispherical ends shown in Fig. (5-33a) . View 
factors F ds / F R g/ F ET and the corresponding components of the 

incident thermal energy per orbit are given in Table 5-1. The 
first term in each block of Table 5-la corresponds to the 
lateral surface of the cylinder, while the second term is 


(a) 


(b) 



Cylinder + Hem> sphericet Ends 




Sun-Ortented Solar Panels 


Figure (5 - 33) - Flat, Hemispherical and Cylindrical Spacecraft Surface Configurations. 


*See Footnote **, p.l. 
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TABLE 5 - 1 

(a) View Factors and Normalized Incident Thermal Energy Per Orbit 
For a Cylinder/Hemispherical-End Spacecraft Configuration 


r DS 


-^es <T ? <7 ?es 


W DS < 2r ?eso) 


LV 


LH 


POP 


Definitions 


f r< 

1 

> 

o<tf $ oo°-e 0 

0 s = 90° 

Ice c *s+AIcs C ^s 

I sc + XIss 

Ice C «s + Alcs s 

2 Js 90 + ^I$S 

Jcc^s+^cs ce s 

Iso + AI$s 

WRst2neso/ 3 > 

0 < 0 < 90° 

13 = 90° 

2r cc c(3+2,(2l* c (3) 

2*1^ +2. (2* I*) 

2J CC C P + A (21* °P) 

2 tT Jjgo + 2^( 2;r Is$) 

2J cc cp + X(2I* c(3) 

2 tt Jjq ( 2*1^) 


ET 


-180° <180° 



(•*!?•) 


Coefficient 


F* 



L t i»>/ t ^ ^cc ^ A '"-cs r* x sc H *$$ ' I ■ co ■ ^ s' 


2E(k, T) es )+X(2rj ej ) 2J c( ,c(3 + X (21^ °/3) 2ir J* 0 + X (2*1*) 2ir F * 0 +X (2irF s 


2ij w 00 + x (2 ij*) 


0 S = cot (c/3 raj) , 0 < « s < 180° 
r = C0 Q = R/(R + H), R = earth radius; H = local altitude 
k= cj3 

neso = 90° + e Q . n es = 180° - cos 1 (S V c/j) ; 0= ifof -90° 

(b) View Factor Coefficients ** 


H Arbitrary 


2 <1 -Vi - r 2 ) 


[ sin - 1 r Wl — r 2 /(1 + r + r 2 /2) - (2/*)Vl - r 2 E(r)1 


(2/3r) [(2+ r 3 ) — (2 + r 2 )Vl — r 2 ] 


(2/3rir) ( (2 + r 2 )(1 - r) [K(k r ) - F (tf r , k r )] -(2-^(1 + r) [E(k r ) - E (tf r , k,)] 
- 2r (1 — r 2 ) } ; k r = 4r/{1 + r) 2 , <f> r " sin -1 (Vl + r/2) 


irr 3 /8 + (3/4c) [ sin “ 1 r — r (1 + 2r 2 /3)Vl - r 2 ] 


(1/4jt) { (1/2r)(3+ r 2 ){1 — r 2 ) In [(1 + r)/(1 -r)] - (3 - r 2 - 2r 3 ) } 


(r/4) II - r (1 -Vi - r 2 )] + [(1 + r 2 )/2r] sin _1 r + (Vl-r 2 /4r) [ 2.5 (1 + r + r 2 ) 
+ (r/2) - (1/jr)(4r 2 + 7) E (r) + (2/ir)(1 - r 2 ) K (r)] 


H= 230 NM 


0.326 


g 


1.113 


1.295 




♦The expressions given for F R are exact for the range of 6 5 stated. Except for values of 7 ? and 0 approaching 90° 
however, the expressions for 0°«9 s <90° - 0 o also approximate well the exact case when 90° - 0 O < 0 s < 90°. 
When 0 S > 90°, F R is negligibly small. For typical low earth orbits f 1 80 < H < 240 NM , note also that 18.13° 


< 0 0 < 20 80°. 


•♦Complete Elliptic Integrals are designated by. K(k), K(r) (1st kind); E(k), E{k r ) or E(r) (2nd kind) 
Incomplete Elliptic Integrals are designated by* F(0 r , k r ) (1st kind);E(© , k) or E{<j> r k r ) (2nd kind) 
h *Data obtained numerically, no analytical expression available 































































BELLCOMM, INC. 


77 


associated with the hemispherical ends.* Since the cylinder 
lateral cross-sectional area (PL) is used as the reference area 
(A) , the factor, A_ = ttD/ 4L, accounts for the maximum projected 
area ratio between the ends and the lateral surface . In Table 
5-lb coefficients of the view factors in Table 5-la are tabu- 
lated: (a) as functions of the altitude parameter, r=R/(R+H), 

and (b) for a specific altitude, H=230 NM. 

The variation in incident thermal radiation 
(F S =Q S /S A ** and F ET =Q ET /SA) is shown in Fig. (5-3 4a) as a 

function of ft and orbital position n for the POP mode at an 
orbital altitude of 230 NM and D/L=0.2. The relationship 
between incident radiation in LV, LH and POP attitude modes is 
shown in Figs. (5-34b) and (5-34c) for ft=0° and 60®. The^ 
corresponding incident thermal energy per orbit (W g and W^) 

is shown in Fig. (5- 35a) as a function of ft for all three modes. 
These results show that the incident energy obtained in POP 
exceeds that in LV and LH for |ft| s 36°. However, the incident 
energy obtained in LV and LH increases markedly as the earth 
shadow interval disappears (ft>69.6°). 

The effect of the hemispherical-end surfaces can be 
assessed by comparing the incident thermal energy data with 
that shown in Fig. (5-35b) for a cylinder/f lat-end configuration 
See Fig. (5-33b).*** Although the incident thermal energy level 
is somewhat higher due to the larger surface area of the hemi- 
spherical ends , the basic trend is the same for all three atti- 
tude modes. Results for other end surface configurations lying 
inside the hemispheric volume should fall within the data for 
hemispherical and flat-end cylinders. 


*In this example only the total incident thermal energy on 
the configuration is considered. Consequently, the hemispheri- 
cal ends are equivalent to a single sphere of the same diameter 
insofar as evaluation of view factors is concerned. In other 
analyses, because of different surface coatings for heat load 
modification, the incident radiation on elemental areas of the 
spacecraft may be of interest. This can be achieved by approxi 
mating the actual surface with an appropriate number of flat 
plates for which the view factors and incident thermal energy 

21 

can be evaluated separately. 

**See Footnote ** on page 7 4 . 

***View factors F Dg , F Rg and F ET for this configuration are 
given in Ref. 21. 



Figure (5 - 34} - Normalized Incident Thermal Radiation vs Location in Orbit { 77 ) 
{LV, LH and POP Modes). 














(a)Spacecraft* Cylinder/Hemispherical Ends (b) Spacecraft Cyfirtder/Flat Ends 




vo 

I 


Figure (5 - 35) - Normalized Incident Thermal Energy Per Orbit vs 0 
(LV, LH and POP Modes) 
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5.4.2 Flat Surface Configuration in Solar Oriented Attitude Mode 

An attitude mode and flat surface configuration of 
interest is one with, a sun-oriented solar panel as encountered, 
for example, in the solar-inertial mode and in the POP mode 
shown in Fig. (5-33c) . Here direct solar radiation impinges only- 
on one side whereas reflected solar and earth thermal radiation 
may impinge on both sides, but varies with 8 and n. The normal- 
ized incident radiation for the sun-oriented panel is shown in 
Fig. (5-36) as a function of 8 and n-* The incident thermal 
energy per orbit vs. 8 is given in Fig. (5-37) . 


*View factors F Dg , F RS and F ET for this configuration are 
given in Ref. 21. 
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Noon 


Location in Orbit, t? (deg) Midnight 


Figure (5 - 36) - Normalized Incident Thermal Radiation 

vs Location in Orbit {Sun-Oriented Solar Panel) 



/? {deg) 


Figure (5 - 37) - Normalized Incident Thermal Energy per Orbit 
vs (3 {Sun-Oriented Sojar Panel) 
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6.0 AVERAGE' EFFECTS OF THE g VARIATION 


For a particular mission the g profile in circular 
orbits is influenced by five parameters: time (t) , time of 

day of launch, (t^) , launch, date (D.^)*, orbit inclination (i) 

and orbit altitude (H) . In previous sections a number of 
system performance factors have been described and related to g 
From a mission design standpoint it is of interest to examine 
mission averages of these factors**- and the effect of mission 
duration (T^) , launch time uncertainty (t L ,D L ) and orbit inclin 

ation and altitude. Before discussing averages of specific 
performance factors, mission averages are defined in general 
and several properties are described. 

6. 1 Mission Averages 


The mission average of a function f(g) can be 
expressed as 

f = f (T M ,t T ,D t , i,H) 
av av M' L' L 


_T 


M 


_ 1 


T 

r m 


T 


M 


f (g)dt = =±- 
'0 M ^ 


f (t,t L ,D L ,i,H) dt 


(6-1) 


0 


where T^ is the mission duration. The extremes on f &v with 
respect to launch conditions (t^,D^) are defined as 


M 


<<f = „>> = Max f a „ = Max ^ | f (t,t L ,D L ,i,H)dt (6-2) 


av ^ av 4 .’ ^ T.. 

^L ,D L V d l m 


; 0 


T 

' M 

<f^„> = Min f aw = Min j f (t,t L ,D L ,i,H)dt (6-3) 


av . av . T m 

^l /D l V d l M 


*See Section 4.0 or Appendix A.l for discussion of launch 

time . 

**The relevance of mission averages to mission design is 
discussed in Section 6.2. 



BELLCOMM, INC. 


83 


Frequently launch, time and date are uncertain at early stages 
of mission planning. If and are regarded as independent, 

uniformly distributed random variables , the mean of f is 

given by 


av 


= Mean 

V D L 


'av 


n r D LM r t LM 

D t I / f av( T M ,t L ,D : 

LM LM J Q J 0 


L ,i,H}dt LM d D L M 


(6-4) 


where 


and 


t T „ = 24 hours 
LM 


(6-5) 


= 365 days (6-6) 

A particular function of interest is f(8) = |b|. In 
Fig. (6-la) the mean and limits of |g| obtained numerically 

by scanning t L and are shown for a 35° , 230 NM circular 

orbit and mission durations up to one year. The mission aver- 
age of | 3 | for all possible launch times and dates lies within 

the envelope formed by <<|b| >> and <|b| >. The [ 3 | „ envel- 

av a v a v 

ope converges to | 3 | =23.53° which is independent of mission dura- 

av 

tion and altitude.* Oscillatory behavior of the envelope is due to 
the multiperiodic nature of 3. The rate of envelope convergence to 

| 3 | is a function of altitude, although the effect is only 
av 

slight in low altitude orbits. Hence, mission average charac- 
teristics for | 3 | , as in Fig. (6-1) will apply in near earth 

orbits even with altitude decay. 


*In calculating | 3 | all possible 3 for a given orbit 
inclination are weighted equally, since t L and are assumed 
to be uniformly distributed. Consequently only orbit inclina- 
tion, which determines the maximum range of 3, has any effect 

on | 0 | . (See section 4.1 for discussion on the 3 envelope.) 

av 
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The variation of the mean and limits of I g I with 

1 1 av 

orbit inclination is shown in Fig. (6-lb)* for a 28 day mission. 

As mission duration increases, the | fi | envelope converges to 

| B | , although the rate of convergence decreases continuously 

as i approaches i , the inclination for a sun-synchronous orbit. 

(See Section 4.4). At i=i the |g| envelope does not converge 

to |si_„/ since $ no longer oscillates between the upper and 

civ 

lower boundaries of the 0 envelope, as shown in Fig. (4-10) . 


*The curves are plotted for 0°<i<180°, although only certain 
inclinations are achievable with a direct launch from a given 
latitude and launch corridor (azimuth constraint) without yaw 
steering during launch or subsequent plane change maneuvers. 

See Eq. (2-9) . 
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6.2 Mission Average of System Performance Factors 


Mission averages of system performance factors are often 
relevant to mission planning, since system designs based solely 
on worst case 8 conditions, which occur only at discrete times, 
may impose excessive requirements or may not utilize available 
performance margins. The system performance factors evaluated in 
Section 5 included: 


1. | ah) /ah - the normalized bias gravity 

r ' [dA gradient torque-impulse on space 
craft in a solar-inertial mode 
(Section 5.1), 


2 . 


‘SL 


the fraction of an orbit in sun- 
light (Section 5.2) , 


3. 


E * 
max 


4. W ** 
s 


the normalized maximum electrical 
energy available from a planar 
solar array, optimally articulated 
(SDOF) with 8, in LV and LH opti- 
mum roll modes, LV and LH cyclic 
roll modes and LV, LH and POP 
fixed roll modes (Section 5.3), 

the normalized direct and 
reflected solar component of 
incident thermal energy per orbit 
on a cylindrical/hemispherical- 
end spacecraft configuration in 
LV, LH and POP attitude modes 
(Section 5.4). 


*The normalization inherent in E „ represents the maximum 

max 

array energy output (E ) available m a particular attitude 
~~ mdx 

mode relative to the output from a sun-oriented array (E ) . 

b 

Since E varies with 8, it is more appropriate in evaluating 

o 

mission averages of the array energy output to normalize rela- 
tive to the sun-oriented array output at some particular 8, say 
8=0°. In view of Footnote *, p.44, 

A A A 

E - E /E = E E 

o max max so s max 

A 

Curves of E_ vs 8 are shown in Fig. (5-32) . 

' o max ^ 

A 

** W ET' t * ie earth thermal component, is not considered here 
since it is constant for all 8 in these attitude modes. 
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The results obtained from evaluating the mean and limits of 
mission averages of these functions, f(3), are presented in 
Figs. (6-2) through. (.6—71 • These results are based on a 230 NM 
circular orbit, although, altitude changes of *30 NM or more 
(in many cases) have only slight effect. The various averages 
are plotted versus mission duration for i=35° and versus orbit 
inclination* for a 28 day mission. As was the case for |p| av / 

the envelope on the functions (f } converges to the mean 
(f ) , which is independent of mission duration, although not 

clV 

necessarily independent of orbit altitude.** 

Mission averages of system performance factors are 
relevant to mission design in assessing certain systems 
requirements and/or performance margins , as will now be 
discussed for the performance factors described above. 

I AH I /AH 
1 — 1 ' max 


For bias momentum dumping via reaction thrust, 
|s2(3+$) | is directly proportional to the propellant 
requirement , 


W = W (T../T) 
p pm M' 


s2 (g+$) 


av 


where W is the maximum propellant consumption per orbit,*** 
T is the orbital period and T^ is the mission duration. The 


*See Footnote *, p. 84. 

**Generally f is a function of altitude, if f(3) depends 
on altitude (H) explicitly, e.g., f(B)=ri es depends on H explicitly 
whereas f (3) = |s2 (3+$) | does not. 

***Since the maximum reaction thrust impulse requirement per 
orbit for dumping | AH | is AH , 

W = (2AH )/I L = T T/I L 
pm max' ' sp gmx ' sp 

where I is the propellant specific impulse, L is the distance 

sp 

between a pair of thrusters fired as a couple and AH is 

max 

defined in Eq. (5-3) . 
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limits on |s2(.8-H) | „ shown in Fig. (6-2) indicate the possible 

variation in the propellant requirement due to launch time 
uncertainty for short duration missions (e.g. T M <90 days) and 

in early phases of long duration missions. Angular displace- 
ment ($>=23°) of the spacecraft principal and geometric axes 

has relatively little effect on |s2(8+$)| aV vs i except at low 

or high inclinations where js2(g+$) | is up to 40% more for 

$=23° than for (j>=0°. 



For solar observation experiments (T ot ) ___ is of inter- 

oij av 

est , since it is proportional to the total-solar-viewing-time 

(TSVT) 3 ' 13 for long duration missions <T M >90 days, say). In 
fact* 

TSVT " <®- 8) 

The limits on (T ot ) „ shown in Fig. (6-3) indicate the possible 

bJj av 

variation in TSVT due to launch time uncertainty for short 
missions. 

A 

E 

o max 


Solar array sizing is usually based on a worst-case 
analysis to insure that basic system power requirements are 
always accommodated. In the absence of constraints a designer 
may be able to afford the luxury of adding all subsystems to 
the basic system requirement without regard to intermittent 
operating requirements of many accessory systems (e.g. experi- 
ments, centrifuge operations, etc.). Since the worst-case 
condition corresponding to certain & angles exists only at 
discrete times however, some margin is available at other times 
to accommodate accessory systems on a scheduling basis. From 
a power availability standpoint the best on-time for accessory 
systems could be selected during a mission by an on-board 
scheduler monitoring the 8 profile and electrical power margins - 
or prior to launch by time-line analysis, when system design 


*The viewing time for experiments affected by earth's 
atmosphere is less. Atmospheric effects on solar viewing time 
are considered in References 13 and 22. 



|S2(jJf23 0 )| 
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and launch time are firm. In the interest of utilizing this 
additional capacity a question for mission designers is: What 

is the mean margin available during a mission for accommodating 
accessory system requirements and the effect on array size? If 
this margin exceeds these requirements, on-time scheduling of 
all ‘accessory systems is a possible design option. If not, 
certain accessory system requirements must be added to basic 
system requirements and array size correspondingly increased. 

The curves in Fig. (5-32) indicate that the worst-case 

condition for all modes considered is 6=0°, where E „ is a 

o max 

minimum. To obtain the same energy output at 3=0° m a parti- 
cular mode, the required array size relative to the sun-oriented 
array size is given by the factor* 


= (E ) . /(E ) . = 1/(E ) . 

s' mm' o max min o max mm 


(6-9) 


Thus, in Figs. (6-4) and (6-5) the corresponding mean and limits 

of (E ) can be scaled by p so that the mean energy margin 

_ o max av J K J 

(E ) in each mode becomes 
mar 


mar 


2 p[(E 


o max av 


) - (E 

' "ITT ' r 


max^min^ 


= p (E ) -1.0 

^ o max av 


( 6 - 10 ) 


This expresses the margin available for accessory systems as a 
fraction of the array energy output at 3=0°^ Similarly for the 
sun-oriented array the mean energy margin (E ) , is 


(E ) = (E ) - (E ) . = (E ) - 1.0 

s ' mar s av 's' mm s'av 


(. 6 - 11 ) 


Mean energy margins and array size factors for the various 
modes are given in Table (6-1) for 35° and 50°, 230 NM orbits. 


*As defined in Footnote *,p.44: E = n /n 

' v s , es / 'eso 
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0 60 120 180 240 


T M (days) 



Figure (6 - 4) - Mean and Limits of (E 0 max ) av vs, Mission Duration 

and Orbit Inclination (LV, LH & POP Fixed Roll Modes) 



0 30 60 90 120 150 180 


t fcJeg) 


A A 

Figure (6 - 5) - (E 0 max ) av and (E s ) av vs Orbit Inclination 
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Table 6-1 Mean Electrical Energy Margins and Array Size Factors 
in 350 and 50°, 230 NM Orbits 


Mode 

A 

^0 max^min 

E 

L mar 

P 

Pr 



1=350 

0 

0 

10 

it 

LV-Optimum Roll 

0.85 

0.08 

0.15 

1.17 

1.11 

1.08 

LV-Cyclic Roll 

0.76 


0 19 

1.31 

| 

1.18 

LH- Optimum Roll 

0.55 

0.28 


1.81 

■M 

1.33 

LH-Cyclic Roll 

0.46 

0.42 

0.58 

2.16 

1 56 

1.43 

LV, LH & POP 
Fixed Roll 

0.52 

0.22 

0 39 

1.93 

1.63 

1.47 

Sun- 

Oriented 

Array 

A 

^s^min 

(^s^mar 




mm 

i=50° 

1.0 

0.03 

0 08 


While the array size is minimum for the sun-oriented 
array, the energy margin is also the least. Hence, comparison 
of other modes with a sun-oriented mode solely on the basis of 
array size* is misleading since the mean array output over a 
mission is higher in the other modes. To obtain array capacity 
comparable to that of any particular mode, the sun-oriented 
array size should be increased by the factor 

Os = * (v - ®s>mar) ‘S' 12 ’ 

Therefore, the array size for that mode relative to the sun- 
oriented array size is 

P r = p/p s (6-13) 

Values of p r for the various modes are given in Table (6-1) . 


*i.e., array size as determined by array output at 6=0° 
(worst case) . 
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In short duration missions and early phases of long 
missions, launch time has an effect on the energy margins. 

This can be observed in Fig. (6-4a) as the difference between 
the limits of CE q max l av and (E q max ) min . As mission duration 

increases, however, these converge to the mean. The effect of 
orbit inclination on the mean margin is shown in Fig. (6— 4b) . 


Thermal control system design for spacecraft is 
frequently based on a worst-case analysis in which the incident 
thermal energy per orbit is either a minimum (cold case) or a 
maximum (hot case) . This establishes maximum internal heating 
or cooling requirements for maintaining spacecraft modules at 
desired temperature levels . Since the worst-case condition 
corresponding to certain 8 angles occurs only at discrete times , 
some margin on heating or cooling requirements exists at other 
times. Because thermal control requirements are directly related 
to electrical energy requirements, any incident thermal energy 
margins are translatable into electrical energy margins, which 
may assist rn accommodating accessory system requirements on an 
on-time scheduling basis as described previously. 


As an example, the relative thermal energy margins 
obtainable in LV, LH and POP attitude modes will be compared 
on the basis of a cold-case design for the cylinder /hemispheri- 
cal-end spacecraft configuration shown in Fig. (5-33a) . The 

curves in Fig. (5- 35a) indicate that W is a minimum at 8=0° for 

s A 

the LV and LH modes , but in the POP mode the minimum for W 

s 

depends on which varies with orbital inclination.* In 

a Hid a 

Fig. (6-6) (W ) . is plotted along with the mean and limits- of 

a o i ill n 

(W s ) av f° r the LV mode. The mean incident energy margin is 
defined as 


(W ) £ (W ) - (W ) . 

s' mar s'av v s' min 


(6-14) 


Due to the normalization inherent in W , Eq. (6-14) 

s 

the mean energy margin as a fraction of the direct 
ponent of incident thermal energy per orbit on the 


expresses 

solar corn- 
lateral 


*See Section 4.1. 
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surface of a solar-oriented cylinder.* Mean margins for the 

three attitude modes are given in Table (6-2) for 35° and 50 °, 

230 NM orbits.** The maximum margin is obtained in the POP 

mode where (W ) is 29% and 34% of (W ) . for the 35° and 
s mar s mm 


50° orbits. 


* 

Table 6-2 Mean Incident Thermal Energy Margins 
in 30° and 50°, 230 IMM Orbits 

Attitude 

Mode 

(W s ) min 

W s ) av 

<w s ) mar 

i=35° 

i=50° 

i=35° 

i=50° 

i=35° 

i=50° 

LV 

1 26 

1.26 

1.35 

1.42 

0.09 

0.16 

LH 

1.16 

1.15 

1.28 

1.36 

0.13 

0 21 

POP 

1.15 

1.06 

1.48 

1.42 

0.33 

0.36 


{Cylinder/Hemispherical— End Spacecraft: D/L=0 2) 


The effect of launch time on incident thermal energy margins m 

short duration missions or early phases of long missions is 

observable in Fig. (6-6) by the difference between the limits of 

(W ) and (W ) . . As mission duration increases, this differ- 

s 9.v s mi n 

ence converges to (W ) . The effect of orbit inclination on 

^ s mar 

A 

(W ) is shown in Fig. (6-7) . 

S clV 

6 . 3 Approximate Evaluation of Mission Averages of System 
Performance Factors 


Because of the number of parameters involved, evalua- 
tion of mission averages (f av .) and subsequent determination of 

the mean and limits is generally a lengthy, although straight- 
forward, numerical process. In the course of this work it was 


*When solar-oriented, the cylinder axis is assumed normal 
to the sun line. 

**Since the curves of W g vs 8 in Fig. (5-35) for the cylmd- 

er/flat-end and cylinder /hemispherical-end spacecraft configur- 
ations are similar, results for (W ) will also be comparable. 
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Mean and Limits of (W s ) av vs Mission Duration 
and Orbit inclination (LV Mode) 
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observed that the mean and limits of f could also be approxi- 
mated with good accuracy in most cases by evaluating f(B) with 
the mean and limits of | B | , that is , 

a V 


f 

a. 

f ( ! 6 

Uv* = f B 

(6-15) 

3V 


<f > 

'V 

f (< 1 

Slav'* 3 'V 

(6-16) 

av 

A/ 

< <f > > 
av 

A/ 

f (<< 

l*lav >> > 3 <<f 6” 

(6-17) 


For comparison, mission averages of several performance factors 
evaluated both ways, are tabulated in Table 6-3. 


Table 6-3 Comparison of Exact and Approximate Evaluation of 
Mission Averages of Typical Performance Factors 

Inclination 

A 

3 

V 

l*W 

«!j8|3V» 

f(/3) 

A 

or* 

< 

V 

Exact 

«v> 

A 

i*- 

V 

Ipproximat 

H 

5 

« f „» 





MPMI 

0 83 

0 87 

0 91 

0 83 

0 89 

0 92 






0 52 

0 63 

0 77 

0 51 

0 60 

0 75 

i=35° 

10 400 

23 56° 

38 07° 


1 28 



■9 

1 34 






Tsi 

0 62 



| 

0 62 






|s2^| 

0 34 

0 62 

0 84 

0 35 

0 73 

0 97 






0 85 

0 90 

0 99 

0 85 

0 91 

0 94 






0 55 

0 72 

0 97 

0 54 

0 67 

0 88 

i=50° 

15 01° 

30 77° 

48 63° 

A 

w s *»* 


142 

1 62 

1 30 

1 38 

1 51 





T S l 

0 62 

0 66 

0 74 

0 62 

0 63 

0 68 





!s2/3| 

048 

0 69 

0 85 

„ 0 50 

0 88 

0 99 


•LV-Cychc Roll Mode **LV,LH & POP Fixed Roll Modes **«LV Mode 
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The approximate method was found to be reasonably 
accurate when f ( B ) was a weakly nonlinear function of g. The 
approximation is less accurate when f(g) is highly nonlinear 

in the range of £ encountered (e.g., f ( B) = [ s 2 B | where |g| >45°) 

Since I 8 I r has been obtained once and for all as a function of 
1 1 av 

orbit inclination* in Fig. (6-lb) , approximate evaluation of f 

a v 

for other orbit altitudes or for performance factors other than 
those studied here, is simple, once f(g) is obtained, either 
analytically or numerically. Convergence of the limits on |&l av 

to I 3 1 is altitude dependent in general, although not signifi- 
av 

cantly different from the curves in Fig. (6-1) for altitudes in 
the range, 230 *30 NM. Hence, approximate evaluation of the 
limits on f &v based on these curves is also reasonable for alti- 
tudes in this range. 


av 


*As noted previously, |g 


is independent of orbit altitude 
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7.0 SUMMARY AMD CONCLUSIONS 

Orientation of the solar vector relative to an earth 
orbit has been defined in terms of two pointing angles : £ , 

the minimum angle between the solar vector and the orbital plane, 
and ip, the location of orbital noon relative to the ascending 
node. Analytical expressions for evaluating 3 and ip have been 
derived as a function of four angles, two that are constant: 

i - the orbit inclination relative to the equatorial 
plane* , 

e - the angle between the equatorial and ecliptic 
planes (23.45°), 

and two that vary with time: 

Y - which describes the earths' motion about the sun, and 

G r - which describes motion of the orbit ascending node in 

1 the equatorial plane due to earth oblateness effects. 

In general, 6 varies harmonically within an envelope which depends 
on both y and i as shown in Pig. (4-3) for both posigrade (0°<i<90°) 
and retrograde (90°<i<L80°) orbits. The angle ip increases mon- 
otonically with time except for orbit inclinations (90°-e<i<90°+e) 
where the orbital plane may pass through an orientation normal 
to the sunline and the monotonic variation occurs only at inter- 
vals . Various- characteristics of the solar pointing angles are 
described in Section 4 and typical time variations are shown m 
Figs. (4-1), (4-2, (4-4) and (4-8). A general summary of properties 
is given in Section 4.5 together with specific results for 35° 
and 50°, 230 NM circular orbits, which have been considered for 
AAP Workshop missions. 

The relationship between a number of system performance 
factors and the solar vector orientation relative to the orbital 
plane has been evaluated for circular orbits. Specific perfor- 
mance factors considered were: 

1) the bias gravity-gradient torque- impulse on space- 
craft m a solar- inertial mode (Section 5.1) , 

2) the orbit sunlight interval (Section 5.2) , 


*Periodic variation of orbit inclination due to earth 

oblateness is a relatively small effect 6 and is neglected in 
this work. 
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3) the electrical energy per orbit available 
from solar arrays in LV, LH and POP modes 
with incrementally articulated* solar 
arrays and various spacecraft roll profiles 
(Section 5.3), 

4) the incident thermal energy per orbit on a 
cylindrical spacecraft in LV, LH and POP 
modes. (Section 5.4) 

While possible options are virtually endless, spacecraft and 
solar array models and these attitude modes have been chosen 
to illustrate typical properties and effects of the B variation, 
as they may effect AAP and future space station missions. All 
results are given in normalized form to facilitate application 
or approximation to particular spacecraft configurations or 
changes in orbit altitude where appropriate**. In addition to 
the analytical formulation these relationships are plotted vs 
5 for 230 NM orbits and other parameters relative to each case. 

Included in Section 5 . 3 are curves showing the maximum 
array energy output relative to a sun-oriented array for each 
mode with optimum array articulation as a function of £. The 
highest array energy output is obtained in a LV-Optimum Roll 
mode where the output is always within 15% of the sun-oriented 
array output (worst case: 3=0°). Because the spacecraft is 
rolled such that the array articulation axis remains normal to 
the plane formed by the sun line and local vertical, large roll 
maneuvers (up to 180°) must be executed abruptly near 3=0°. An 
easily implemented, suboptimal alternative at low 3 is the LV- 
Cyclic Roll mode in which array output is within 23% of the 
sun-oriented array output up to 3=60°. Array performance with 
comparable roll profiles in LH modes is decidedly lower than in 
LV modes for |&|<45°. A more detailed summary and comparison 
of optimum array performance in various modes is given in 
Section 5.3.5. 

Mission averages of system performance factors are 
often relevant to mission planning, since system designs based 
solely on worst case 3 conditions, which occur only at discrete 


*Solar arrays are assumed to be articulated about one axis 
and only intermittently (e.g., daily), not at rates comparable 
to the orbital angular velocity. 

**Although the analysis is based on circular orbits, the 
results are also applicable to slightly elliptical orbits of 
comparable period. (See Footnote *, p.4). The equivalent circu- 
lar altitude (H) is the average of apogee and perigee altitudes, 

H= (h +h )/2. For arbitrary elliptical orbits, however, system 
a p 

performance factors are a function of both 3 and tj>. 
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times, may impose excessive requirements or may not utilize 
available performance margins. Mission averages (f av ) * of 

each performance factor, (1-4) above, have been evaluated as a 
means of assessing the effect of the B variation and arbitrary 
launch time and date. The mean and limits of f based on 
independent, uniformly distributed launch time and date are 
plotted vs mission duration and orbit inclination in Figs. (6-2) 
through (6-7) . The limits of f . which asymptotically 

a. v 

approach the mean, (f ) , as mission duration increases, indicate 
the possible variation in f due to launch time uncertainty. 

In evaluating whether or not a design based on worst 
case conditions (B ) is entirely realistic, the difference 

(or margin) between f (B ) and the mean or limits of f &v can be 

useful data. In one case studied in Section 6.2, f ( 3 > = | s2 3 1 

is proportional to the propellant requirement per orbit for 

momentum dumping of the bias gravity gradient tor que- impulse by 

reaction thrust. Here, f(f3 )=f(45°)=1.0 and f =0.62 for a 35° 

' wc av 

orbit. Hence propellant requirements based on the worst case 

over a long mission (e.g. , 56 days) would be excessive. If 

requirements are bas ed n evertheless on the worst case, the 

calculation, f(B )-f . indicates a propellant margin up to 38%, 

wc av 

which might be used for supporting short duration experiment 
pointing modes. In Section 6.2 mission averages of the perform- 
ance factors, (3-4) above, are used to evaluate the margins on 
electrical energy per orbit from solar arrays and the incident 
thermal energy per orbit on spacecraft in various attitude modes. 
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*The mission average (f __) of a function f(g) 

" av 

as the time average of f (6) over a mission. 


is defined 
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APPENDIX A 


Location of Sun Line and Ascending Node at Launch 


The purpose of this appendix is to determine as a 
function of launch time the initial location of the sun line 
and the orbit-equatorial plane ascending node as represented 
by the angles y and fl ro . Since launch time depends in 

general on the reference point for measuring time, some 
discussion with regard to this point is in order. In 
Fig. (A-l) a particular orientation of the ecliptic, equator- 
ial and orbital planes is shown, similar to that shown in 
Fig. (2-2) . In the equatorial plane, true noon is defined by 
the intersection with the equator of the meridian plane contain- 
ing the sun line and the normal to the equatorial plane 
(earth's polar axis). With an analogous definition of noon 



Figure (A • 1) - Orientation of the Orbital Plane Relative to the Equatorial and Ecliptic Planes at Launch 





bellcomm, Inc. 

Appendix A 


' A- 2 


for the orbital and ecliptic planes it follows, because of the 

relative orbit inclinations, that true noon occurs at different 

points for the three planes as noted by N , N and N in 

2 o eg ec 

Fig. (A-l) . During orbital flight; the logical reference would 
be N q . For launch purposes, however, N is appropriate. 

A. 1 Location of Sun Line 


In Fig. (A- 2) launch time is measured from two 
reference points , one based on autumnal equinox and the other 
ll 

on midnight (00 ; 00 ) of September 23. The various time intervals 
are related by 

t + 24d T = 24 D t + t T (A-l) 

510 -b J_i 

so that 

= D l + ‘VW 724 (a - 2) 

where 

d L = launch time measured from autumnal equinox in 
days (non-integer) 

D = number of calendar days (integer) at launch 
since September 23 (00:00 U.T.*) 

t T = launch site Universal Time* at launch 
(0<t L <24.0 hours) 

t = launch site Universal Time* at autumnal equinox 
ae (0<t <24.0 hours)** 

d6 


*Here Universal Time (U.T.) is defined as Greenwich Mean 
Time corrected for the longitudinal location of the launch site. 

For KSC, located at 80.5° W. Longitude, the launch site correction 
is -80.5/15 = -5.37 hours. 

23 

**Generally t is obtained from emphemeris data. However, 

ac 

the maximum error incurred in y by ignoring t T and t is less 

, O JLi 3.0 

than 1°. 
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Upon evaluating d_ , the sun line location at launch (y ) can be 

L 23 ° 

determined exactly from ephemeris tables or approximately from* 


y 


o 


(A- 3) 


where y is the mean sun 


line rotation rate as given m Eq. (2-8) . 



24D l 


l 

x ae 

24<I l ' 



I > 1 


00 h 00 J 
Sept 23 

1 * < 

k 00* 1 00 00 h 00 J 
Sept 24 Launch 

Day 

j - - 

1 00 h 00 
Launch 
Day + 1 


Universal Time 
at Launch Site 
(Hours) 


Autumnal Launch 

Equinox 


Figure {A - 2) - Relationship of Launch Time to Autumnal Equinox 


A. 2 Location of Ascending Node 

In Fig. (A-l) , the point L is considered to be 
immediately above the launch site. The longitude of the 
launch site with respect to the noon meridian (through N ) 
can be specified by the angle p T where g 

i-i 

P L = [15(t L +t N ) - 180°] (A- 4) 


t L is launch site Universal Time as defined in Section A.l, and 

t N is the difference between the true noon and apparent (clock) 

noon at Greenwich.** As is evident from Fig. (A-l) , the longitude 
of the launch site can also be specified with respect to the 
autumnal equinox reference (Z ) by either of the angles. 


*See Section 2.2 and Footnote *, p. 7. 

**True noon represents (t L +t N )-12.0. The variation in t^ 

23 

is available from ephemeris data. Since |t N ( max <0.27 hours, 
the maximum error incurred in p T (and Q ) by ignoring t„ is 4°. 

JLi 3TO IN 
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(^ ro +ft L ) or . Consequently, the initial location of the 

ascending node is given by 


£2 ~ Pt + Jin — £2 T 

ro L N L 


(A- 5) 


where £2 N and & L are defined in Fig. (A-l) and derived m 

Sections A. 3 and A. 4 respectively. Substituting Eqs. (A-4) , 
(A-9) and (A- 15) into Eq. (A-5) yields 


£2 = [15 (t T +t„) -180° ] + tan 1 (tany ce) - tan ^ (tana T s£) (A-6) 

ro L N ' 'o L 


where y Q is the sun line location at launch, e = 23.45°, is 
the launch azimuth and £. is the latitude of the launch site. 


A. 2.1 Determination of £2^ 

As shown in Fig. (A-l) , £2 N is the' angle between the 
noon meridian plane (through N ) and the autumnal equinox 
reference (Z^). Consider a coordinate system (x^ ^ defined 
such that y N is parallel to Y £ (the normal to the equatorial 
plane) , z^ lies in the equatorial plane along the line of inter- 
section with the noon meridian plane and Xj^ completes the right 

hand system. The transformation relating (x N ,y N ,z N ) and 

(X , Y , Z ) for a given launch time (y ) is: 

S S S O 




= a £s 






(A- 7) 


Since x^ is normal to the noon meridian plane, it follows that 
the element a^g, which is the direction cosine between x^j and 
Z s (the sun line) must be zero. This yields 



so that 


- c£2. 7 sy ce - s£2 M cy^ = 0 
N 1 o N o 


(A- 8) 


tan£2 M = tanY • ce 
N o 


(A-9) 
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A. 2. 2 Determination of Launch Parameters {S2 T , i) 

Jj 

The relationship between various launch parameters 
can be described in terms of the diagram in Fig.(A-3). 


Earth's Polar Axis 



a l Launch Azimuth 
^ — Launch Site Latitude 


i — Orbit Inclination 

12 L — Launch Site Equatorial Longitude 
(Relative to Ascending Node) 

12 0 — Launch Site Orbit Longitude 
(Relative to Ascending Node) 


Figure (A - 3) - Launch Site Coordinates and Parameters 

The coordinate system ( x n 'Y n r z n ) is the nodal system 
used previously. The coordinate system (x T ,y T ,z T ) is defined 
such that z T corresponds to the launch site local vertical, x T 

-b -L 

is parallel to the launch heading and y L completes the right 

hand system and is normal to the orbital plane. 

The launch azimuth a L is defined positive with 

respect to geographical north with a T = 90° corresponding to a 

due east launch. The transformation relating {x L ,y L ,z L ) and 

(x n ,y n ,z n ) coordinates can be expressed in either of two 

elementary Euler angle sequences: 90°-a T ) or (i, £2 ) 

so that L L o 



B — n 


= Cy 
— n 


(A-10) 
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where 


and 


B = 


m ^ T Y 

90 °~a L l -z 


(A— 11) 


(sa^cfij^-cctj^sJlsfl^) 

Ca L CJl 

(-sa^sfl^-ca 

- (cot T c£h +sa r S£sJ2 t ) 

Lt 1j Jj 

sa^cZ 

(ca T -sa. 

J-l Jj ll 

clstt r 

\ *-> 

sZ 

ClCQ.^ 


C = ' 

T y T z 

& i 






V 

r 




(A-12) 


! cfl ci 

! O 

cG si 
o 

-sft 

o 



= 

-si 

ci 

o 




s£2 ci 
o 

sfi si 
o 

cQ 

o 



Of the many possible 

relations obtainable 

from equating ele- 

ments of B and 

C the 

two of 

interest here 

are E ) 22~ c 22 an< ^ 

b 23 =C 23 30 that 







ci = 

= sa^ cZ 

(A— 13) 

and 








0 = ca^sO^-sa L cfi^sii. 

(A-14) 


or 


tan £i L = tan a^sZ 


(A- 15) 
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A. 3 Initial Solar Vector Orientation C 0 Q r ^ 0 ) 

The pointing angles, 6^ and iJ> Q , represent the solar 

vector orientation relative to an earth orbit at launch. These 
may be evaluated in terms of launch parameters (a^, Z, p^) and 

the equatorial sun declination (-a^) at launch. From Eqs. (3-5) 

and (3-6) the expressions for S Q and are 


s$ Q = si ( sy q ce cG ro ~ cy g s ^ rQ ) + ci S Y 0 se (A- 16) 


tan # 0 = ci(sY 0 ce cQ ro - cy q s^ rQ ) " si sy q se 

C Y 0 cfl ro + s Yo ce s \o 


(A-17) 


where 


fi ro ” P L + 


(A-18) 


and (p T , Q t ) are defined in Section A. 2. As given in 

Jj JN lj 

Eq. (4-8) , 


sa, = s y„ se 
d o 


(A- 19) 


so that 


ca. 


I 2 7 2 2 

= *V C Y 0 + S Y„ c e 


d ’o ~ ‘o 

which means, in view of Eq. (A-9) , that 


(A- 2 0 ) 


s& M = sy ce/ca. 
N o d 


(A- 21) 


cSl T = cy /ca, 
N 1 o' d 


(A— 22 ) 
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Substituting these equations 
after some manipulation 


se o = s« d ci 


tan i h = ca ci 
o a 


into Eqs. (A-16) 


+ s ca d 

s ( fl L~ p L ) ~ Sa d 
ca d c ( '** p ) 


and (A-17) yields 

si (A-23) 

si 

— (A- 2 4 ) 


The orbit inclination (i) and the launch site longi- 
tude (fl T ) relative to the orbit ascending node are related to 
Ju 

launch azimuth (a_) and launch site latitude (i) in Eqs. (A-13) 
and (A-15) . The launch site longitude (p L ) relative to true 
noon (N e g) is obtained from Eq. (A- 4) . 

Certain special cases are of interest. For true 
noon (p=0°) and true midnight (p T =180°) launches, it follows, 

Jj 1 j 

in general , that 


= 


sa d ci ± sft L si ca d 


sa^ s ( a d + Si ) 


Noon{+), Midnight (-) (A-25) 


and 


tan \Si = ca-, ci s£2 T + sa^, si 
T o d L - a 

ca d cflj^ 


Noon (+) , Midnight (-) 


s a L cl s{a d +£) - sa d 
Cct d Cot L 


(Noon) 


(A— 26a) 


-s a L cSL s(a^-H) + sa^ 
^ a d Ca L 


(Midnight) 


(A-26b) 


where the latter forms in each case result from substituting for 
ot L from Eq. (A-15) . 
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As indicated by Fig. (A-l) , the angle (n T -p T ) locates 

i-l Lt 

the equatorial noon meridian plane relative to the orbit ascend- 
ing node. Hence, whenever (^ l -p l )=± 90° , it follows that 

sB Q = s(a d + i) (A-27) 

*o * ± 90 ° (A— 28) 

which implies that is at the upper or lower 6 envelope 
boundary.* The corresponding launch time (t T ) at the launch 
site is determined from Eqs. (A-4) and (A-15) which yield 

P L = 15(t L +t N ) - 180° = Q l -(± 90°) 
or 

_(A-29) 


t — — t + 

L TJ 15 


tan ^ (tan a T sZ) + 90‘ 

Xj — 


*See Section 4.1. Algebraic signs (+) are associated with 
upper (+) and lower (-) envelope, boundaries . 
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Gravity Gradient Torque in Solar- Inertial Mode 


In the solar-inertial mode the orxentation of the 
spacecraft geometric axes (x ,y ,z ) relative to the orbital 

SO SC 5 0 

plane is specified by the solar pointing angles (6,ip) shown in 
Fig. (B-l) . The coordinate axes (Xp,Yp f Zp) are principal axes 

where X is parallel to the roll axis (x ) with (Y , Z ) dis- 
P sc D D 


placed by a fixed angle <f> from 

*n • Vlv 



^sc ,z sc* ' 



Satum-V Workshop Configuration 
(ATM Arrays Omitted) 


Figure (B • 1) - Spacecraft Orientation in Sofar Pointing Mode 


The transformation relating 
(x n ,y n ,z n ) is given by 


(X , Y ,Z ) and 
P P P 


p = T X T X T y p = T , X - . T y p 
-p 4. g f — n (3+4.) ip -n 


(B-l) 


* See Footnote ** on p. 34. 
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The gravity gradient torque expressed in general 

form is 12 


= 3n 0 (p,xl.p) (B-2) 

where /y~ = /i/R^~ is the instantaneous orbit angular velocity, 

M is the earth's gravitational parameter, is the local geocentric 

altitude, P is a unit vector directed along the local vertical and 
I is the spacecraft inertia tensor. In Fig. (B-l) the coordinates 
( X £ V /Y£ V / X £ V ^ represent the local vertical system so that 

-Zv ~ The transformation relating < x £v 'Y £v /Z Av ) and 

(x ,y ,z ) coordinates is 
n fj: n n 


—Zv 



in 


(B-3) 


In view of Eqs. (B-l) and (B-3) the unit vector p_ m Eq. (B-2) 
may be expressed in principal axis coordinates as 


p 

r-sc 


= P- 


= T 


x 


ipY T Y p 


3+<j> ^ ~£v 


e+<l> tp-nl 


-s{*-n) 

s (e+$) c ( ^ — n ) 
c (g+$)c (t-n)i 


Substitution of Eq. (B-4) into Eq. (B-2) yields 


T = 3y 


[ ( VV P 2 P 3l 

o \ (I x _I Z )p l p 3 


^ (I y I x )p l p 2 


(B-4) 


(B-5 ) 
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(I -I )s2(g+^)c (ijj-n) 

£. y 

{I -I )c(g+$)s2(+-ri) 

“ (I -I x ) s (3+<}>)s2(iJ;-r|) 


T s2(3+<f>)c (ti-n) 
gmx y 

T c(B+$)s2 <ip- 
gmy r 

T s (3+$)s2(ip- 
gmz ' VT 


2 

For a circular orbit, y = is constant and n = to t where to 

o o o o 

represents the orbital angular velocity. 


< C < £T 
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Evaluation of Electrical Energy Available from a Solar 
Array on Spacecraft in LV, LH and POP Attitude Modes 


The purpose of this appendix is to obtain analytical 
expressions for the electrical energy available from the solar 
array configuration described in Section 5.3 for LV, LH and 
POP spacecraft attitude modes and various roll profiles. 

The profiles include optimum roll, continuous roll, 
cyclic roll, and fixed roll. In each case, it is first 
necessary to determine the array sunlight interval by compar- 
ing the earth shadow boundary (*ri ) with the panel self— 

0S 

shadow boundary (n , ,n ■,) defined in Section (5.3). With 

SS-LSS^ A 

this information, the normalized electrical energy (E)* is 

then evaluated after integration of Eq. (5-16) . 

C.l hV-Optimum Roll Mode 

In the optimum roll case, the panel self— shadow 
boundary is obtained after setting cX LV =0 in Eq. (5—18) with <j> 

defined in Eq. (5-20) by 


tan<j> = -sn/tanS (C-l) 

From this it follows that 

tang 

c<j> « , 5 5 — (C— 2) 

~J tan g + s n 


Rewriting Eq. (5-18) and substituting tan<{> and c<j> above yields 


cX T . T = so cgcri + ca c<j>sg [l+tan<J> (-sn/tang) ] 

i-r V p p 


CC-3) 


= SttpCgcn + C0tpSe/c4. = cgjjsapCn + ca^ *^tan 2 g + s 2 nj 


*Hereafter, the term electrical energy in this appendix 

A 

will imply normalized electrical energy (E) . 
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Equating the bracketed term to zero, transposing terms and 
squaring yields 

s 2 a c 2 n = c 2 a [tan 2 g + s 2 n] 

Jr ir 

= c 2 ot tan 2 g + c 2 a (1-c 2 ti) (C-4) 

P P 

or 

c 2 n = c 2 a (l+tan 2 8) = c 2 a /c 2 8 (C-5) 

P P 

Since cn is an even function, the values of n. obtained from 
Eq. (C-5) are symmetrical about n=0°. Consequently, the self- 
shadow boundary is also symmetrical about n=180°, so that * 


’180° — cos ^(ca /eg) g <_ a 

P 3? 


n ss2 n ssl ^ss 


(C-6) 


.no self-shadow 


8 > a 


A plot of n vs. g obtained from Eq. (C-6) for various values 
of a is shown in Fig. (C-l) . The dashed curve represents the 

ir 

earth shadow boundary (n _) , as obtained from Eqs. (5-6) and 

(5-7) for a 230 NM circular orbit.** The sunlight interval is 
then determined by noting for specific and 8, whether self- 

shadowing occurs first (n <n ) or whether earth shadow occurs 

S S 0 S 

first (n <ri ) . Since the self-shadow boundary is symmetrical 

0 S 5 S 

about n=0 , the sunlight interval (ri-^nir^) is 


* Of the two solutions to Eq. (C-5) , namely cr^ca^/cg and 
cn= -cc«p/c 8, only the latter is of interest, since the other 
solution does not satisfy Eq. (C-3) for cX L ^=0 with a positive 
radical. 

** For orbit altitudes other than 230 NM the curves of 
n. vs. 8 are similar as indicated by Fig. (5-4) . 
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given by 


n 2 




l si 



( n < n 
es s 

no self- 


or 
ss 

shadow 


(C-7) 


if 


n > n 
es 'ss 


180 


(deg.) 


120 



30 60 

0(deg.) 


Legend 

c » eA ■ ^<v m 

cas * B ED 


Figure (C - 1) - Variation in Panel Self-Shadow Boundary with a p and j3 (LV - Optimum Roll Mode). 


iFi view of Fig. (C— 1) , it follows that p „ = n in region B 

s e S 

whereas n s£ =n ss in region A. 

The electrical energy (E) available from the array 
can be evaluated after integration of Eq. (5-16) . Upon substi- 
tuting Eq. (C-3) for cA the integration yields 
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Appendix C 


E = 


es J 


cX LV dn 2n 


_ eg 


2n 


es 


{ 2sa s n 
p 


es 


Ca 

. + 2 — £. 
si eg 


fSJl 

J ■ 

-n, „ 


cX LV d " 


s£ 

‘si 


Vs 2 g + c^gs^n dn V (C-8) 


= — icgsa 
n es \ 


^ n s l 


Sn . + ca 

p s J, p 


~\J 1-c^ gc^n dn^ 


The last term can be recognized as a form of an elliptic inte- 
gral of the second kind, which can be converted to standard 
form by the substitution, n=6-90°, so that * 



where 



gs e de 


E (k , 0+180° ) -E (k) 
E(k,0) + E (k) 


,0+180° 


Vi 


2 2 
-k s 6 d0 


(C-9) 


k=cg 


(C-10) 


and 


0 =n s r 9O< 


(C-ll) 


* The incomplete elliptic integral of the second kind, 

E(k,0), is available in tabulated form^ as a function of the 
modulus k and argument 0. The second function E(k) is the 
complete elliptic integral of the second kind (0=90°). The last 
step in Eq. (C-9) follows from the property that 

E(k,nir+0) = 2nE (k) +E (k , 0) 
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The expression for E can then be written as 


E = -i- (sa 
n es \ 


p c6sn sjl + ca p [E(k,0)+E(k) ] 


) 


(C-12) 


In Fig. (5-9) E vs. 8 is plotted for various panel 
angles (a ) . The upper envelope boundary on the curves repre- 

hr * 

sents the maximum electrical energy (E ) available from the 

max 

array in the LV-Optimum Roll mode with optimum variation in 

A 

panel angle as a function of 8. Plots of E and a , . . 

3 max p(opt) 

obtained numerically are shown in Figs. (5-28) and (5-30). In 

general, an analytical expression for a p ( opt ) can be determined 

by setting 3E/3a =0, although the result may not be amenable to 

P 

closed form solution if n =n 

s £. ss‘ 

In view of Figs. (5-30) and (C-l) however, it follows that 
n s£ =n es ^ or 6 in t ^ ie vicinity of a p =a p ( 0 pt) t ^ ie opti™ 1110 

roll case. Consequently, n in Eq. (C-12) is not a function of 

S Kf 

A 

in that region, so that 3E/3a p =0 is simply 


since n is a function of a . 
ss p 


3E = _JL_ f 
3 “p n es l 


ca cBsn - sa[E (k,0)+E(k) ] ) = 0 
p es p 


I- 


(C-13) 


The optimum panel angle is then given by 

c8sn 


tana 


es 


p(opt) E(k)+E(k,n es -90°j 


(C-14 ) 


which is in agreement with the curve of a p( 0 p t ) determined 
numerically in Fig. (5-30) . Since 


ca 


[E (k) + E(k,n es -90°) ] (C-15) 

V[E(k) + E (k, n es -90 o ) ] 2 + c 2 8s 2 n eg 


P (opt) 
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it follows from Eq. (C-12) that E is given by 

max 3 * 


E = cot 
max 


p(opt) [E(k)+E(k ^es- 90 ° )] * C1 + tan “p(opt) ] ^es 


= [E(k) + E(k,n es -90o)]/n es ca p(opt) 


(C-16) 


n es E(k) + E(k,n_-90°) ] * 2 + c 2 6s 2 n 


es 


es 


which agrees with the results in Fig. (5-28) . 


C . 2 LV-Continuous Roll Mode 

In the continuous roll case, the panel self-shadow 
boundary is obtained after setting ci LV =0 in Eq. (5-18) with 

<j>= -r\. The expression for cA T , 7 with <p= -n is 

J-iV 

C *LV = S ^ a p + ^ cn + ca p cgs 2 n 

2 

= -[cot cgc n ~ s(a +$)cn - ca c&] (C-17) 

P P P 


The solution of Eq. (C-17) for cA^=0, which satisfies the 
constraint, CK|cn|£l/ is given by 


cn 


s(a + B) - -yjs 2 (a +^) + 4c 2 a 


c 2 B 


2ca p cB 


(C-18) 


Since the cosine function is even, the two values of n result- 
ing from Eq. (C-18) , which specify the self-shadow boundary 
(n ccl and n ,) , are symmetric about n=0° sothat t 0 =-n . =n 

bbJ - 5SZ SSZ SSl SS. 

A plot of n vs. 6 obtained from Eq. (C-18) for various values 
of “p is shown in Fig. (C-2) . The dashed curve represents the 
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earth shadow boundary as obtained from Eqs.(5-6) and (5-7) for 
a 230 NM circular orbit. 


180 


150 


(deg.) 


120 



30 60 

0(deg.) 


Legend: 


CaseA 
Case B 





Figure (C - 2) - Variation in Panel Self-Shadow Boundary with a p and (3 (LV - Continuous Roll Mode) . 

Since the self-shadow boundary is symmetrical about 
n=0°, the sunlight interval (-n <n<n e , ) is also determined as 

S )6 S X 

in Eq. (C-7) . In view of Fig. (C-2) it follows that n =n 

s i es 

in region B whereas n „=n in region A . 

S a/ s s 

A 

The electrical energy (E) available from the array 
can be evaluated after integration of Eq. (5-16) . Substitution 
of Eq. (C-17) for cA and integration yields 


-•t/ 


s £ 


2n j_ cX LV dn 
” n s£ 


ca cB 


2n 


i [ 2s (a p +B)sn s«, + 2 (2 

e s 


S £ 


(C-19 ) 
s2n st>] 
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In Fig. (5-11) E vs g is plotted for various panel angles (a^) . 
The upper envelope boundary represents the maximum electrical 

A 

energy (E ) available from the array in the continuous roll 
max ^ 

mode with optimum variation of a with 6. Plots of E and 

p itlciX 

°p(opt) °kt a; i- ne< ^ numerically are shown in Fxgs. (5-28) and (5-30) 


The 'slight bend in the curve of 1 a p ( opt ) vs. g at g=8 T =47 0 signi- 
fies the transition point for calculating E , since n =ri 

^ IUclX 5)o So 

when g<8 . For 8>8 , setting 3E/3a =0 to obtain a , leads 


p (opt) 


leads 


to an unwieldy transcendental equation, since n C!0 =n__ is a 

function of a . The result is not amenable to closed form 
P 

solution. For 8<8 , determination of a , . . is relatively 

t p ( opt ) 

simple since n. cl , = n oc . is not a function of a . As a result 

Id Xj fci fcx fv 


JT = { (2cBsn es )ca p " S“ p 12sesn es + x-C 2 n es -s2n es ) 3 

p es V r c- 


Setting SE/Sa^ = 0 yields 


(C-20) 


tana 


p (opt) tang + K(g) 




(C-21) 


where 


K(g) - < 2 n es -s2n es )/4sn es 


(C-22) 


tang ± K (g) 

*p (opt) j 2 

VI + [tang + K (g) ] 


(C— 23) 


It follows from Eqs. (C-21) , (C-23) and (C-19) with ri £=n eg 

that E is given by 
max 3 x 
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E 


c8 


max 


— J2c8sn sot , . + ca , . [2s8sn + ^—(2n -s2n ) 

2n es ^ es p(opt) p(opt) es 2 es 'es' 


cgsn 


n. 


es 
cBsn 


pc p 

— °“p(opt) Itan6 + K <e>ntan a p(opt) + 1] 


es 


n 


es 


{ 


[tang + K ( 8) ] /ca 


p (opt) 


■} 


c8sn 




es J 1 + I tan 6 + K(8)] 2 


6<8. 


es 


CC-24) 


The results in Eqs. (C-24) and (C-21) for E and a , . . are 

max p Copt) 

consistent with the numerical results in Figs. (5-28) and (5-30) . 


C.3 LV-Cyclic Roll Mode 

In the cyclic roll case, the panel self-shadow bound- 
ary is obtained by setting 0X^=0 in Eq. (5-18) with 4 >=( l ) C y C as 

defined in Eq. (5-21) where 


eye 


+ 180° 

-180° 

< 

T) 

<_ -90° 

-n 

1 

VD 

O 

o 

< 

n 

£ 90° 

- 180° 

O 

o 

< 

n 

< 180° 


(C-25) 


The expressions for cX LV with c are then 


(cx LV ) = s( a +8)cn + ca ces n | n [ <9 0 * 
P p 


and 


(cX TXr ) » s(a -8)cn + ca c8s 2 n 
LV o P P 


= -[ca c8c r\ - s(a ~8)cn ~ ca c8] 
P P P 


(C— 26) 


90° < n I <180° ' (C-27) 
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However, in determining the self -shadow boundary, only the 
latter expression is significant.* The solution of Eq. CC-27) 
for (cX LV ) 2 =0, which satisfies the constraint 0<Jcn|.ll» is 

given by 


cn = 


S( V 6) ~ V s2 <y 


8) + 4c^a 

hr 


2ca eg 
P 


8<a 


(C-28) 


As before the self-shadow boundary is symmetrical about n=0 o 
so that 


180°-cos 


-1 


n ss2 wT) ssl n ss 


-s (g p -8) t-^/s 2 (g p -g)+4c 2 g p c 2 g 


2ca eg 
P 


8<g 
~ P 


no self-shadow 


B> g 

P 

<C-29) 


A plot of n vs. g obtained from Eq. CC-29) for various values 
s s 

of g is shown in Fig. (C-3) . As in previous cases, the dashed 
P 

curve represents the earth shadow boundary for a 230 NM circu- 
lar orbit. 


* Since (cX LV ) ^=0 has no real roots of cn which satisfy 

the constraints 0£[cn|<l and 0°<a , g<90° , it follows that only 

(cX T „) _ need be considered in determining the self-shadow 
XjV a 

boundary. 
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180 
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120 



0 ((tog.) 


Legend 

CaseA H ^ <T ?« 
&M * B □ 


Figure (C - 3) - Variation in Panel Self-Shadow Boundary with a p and 0 (LV - Cyclic Roll Mode) . 


Since the self-shadow boundary is symmetrical about 
n=0? the sunlight interval (-n s ^£n£n s is also determined as 

in Eq. (C-7) . In view of Fig. (C-3) , it follows that n =n 

s a es 

in region B whereas n =n in region A. 

SI s s 

The electrical energy (E) available from the array 
can be evaluated after integration of Eq. (5-16) . Substitution 
of Eqs.(C-26) and (C-27) for cA in the appropriate integration 
intervals yields 
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n 


E = 


2n 


es 


s& 

cx^dn 


l sZ 


( -90' 


2n. 


es 


^■ cX LV^2 dn + 


v - 




90° 

(cX LV } l dT1 
-90° 


f 

y 9Q' 


‘sz \ 

(cA LV ) 2 dn • 

j 


90' 


2n 


es 


(cX LV ) + 2 


f 

90 ' 


s Jl 

2 dn ^ 


ITT" 4ca p S0 + 2s(a p -e)sn sJl + 
es l 


ca cB 

— f— C 2 n sr s2n s 


*’] 


(C-30) 


In Fig. (5-10) E vs. B is plotted for various panel angles ( a p ) ♦ 
The upper envelope boundary represents the maximum electrical 

A 

energy (E ) available from the array in the continuous roll 
max , /. 

mode with optimum variation of a with B. Plots of E and 

p max 

a p (opt) °ktained numerically are shown in Figs. (5-28) and (5-30). 

In view of Figs. (5-30) and (C-3 ) , it is evident that 
n sA =n es ^ or when a p~ a p( 0 pt) *- n t ^ Le c Y c li c roll case. 

Consequently, an analytical solution for a p^ 0 p t j obtained by 
setting 3E/3a =0 with n„ 0 =n. will hold for all B- This 

p S Xj 6S 

procedure yields 


3E 

3a 


(i 


— * (2cBsn „) ca^ 
2n es l es P 


- sa p [4sB-2sBsn es +^-(2n es -s2n es )] j = 0 


(C— 31) 


so that 


tana p (opt) G(6) tanB + K(6) 


(C-3 2) 
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where 


and 


6 ( 6 ) 

K(p) 


2 ~ sn. 


es 


4n 


es 


2t| es - s2n es 


4sn 


es 


CC-33) 


(C-34) 


CCt p (opt) 


G(g)tang + K(g) 

•\/l + fG (3) tang + K(g)] 2 


(C-35) 


It follows from Eqs. (C-32) , (C-35) and (C-30) with r i s ^=n eg 

that E is given by 
max 3 J 


E 


max 


_ / s 0 

J - <2cgsn sa . . v + ca , . J4sg-2sgsn +^( 2 n "S2n ) 


2ri \^ ps,T1 es au p(opt) T ^p(opt) 1 

6S n 


es 
cgsn 


es 1 


es 


n 


es 


ca p(opt) [G(e)tan3 + K(e)][tan a p(opt) + 1] 


C$Sri i 1 2 

— — -yi + [G (g) tang +K (g ) ] 
es * 


n 


(C-36) 


The results in Eqs. CC-36) and (C-32) for E max and “p^ 0 p t j are 

in agreement with the numerical results in Figs. C5-28) and 
(5-30) respectively. 


C.4 Fixed Roll Mode 


(n ssl 
and (f> 


In the fixed roll mode the panel self shadow boundary 
/fL ss2^ '*' s °btained by solving Eq. (5-18) for n with cX^.=0 

^constant. The general form of this equation is 


cX 


LV 


= A + Ben ” Csn = 0 


CC— 37) 
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where 


and 


A = ca c<j> s( 
P 

8 = sap c6 


C — ca eg s<f> 
P 


(C-38) 

(C-39) 

(C-40) 


In contrast to previous cases, the solutions of Eq. (C-37) for 
n are generally not symmetrical about n=0°. The two solutions 
can be expressed in terms of sn and cn after alternately elimi- 
nating cn and sn from Eq. (C-37) and solving the resulting 
quadratic equation for sn and cn respectively. This yields 


and 


where 


sn 


AC + BR 


ss2 


cn 


ss2 


sn 


ssl 


cn 


ssl 


R e -y/? 


“ 8 2 ■ 

t C 2 

... -AB 

+ CR 

B 2 

+ c 2 

_ AC 

- BR 

8 2 

* c 2 

_ -A8 

- CR 

B 2 

+ c 2 

+ C 2 

- A 2 


(C-41) 


(C-42) 


CC-43) 


(C-44) 


/2T 2 27 
/c g - c a c 4> 
P 


(C-45) 



BELLCOMM, INC. 


C-15 


Appendix C 

If 0_<Ca , g)<90°* and — 90°<^<0,** then A and 6 are always posi- 
* ^ n 

trve and C is always negative. In view of Eqs.(C— 42} and 
(C-43) it follows that cn gs2 l0 and Sr i ssl :S 0 so that 90°<_n ss2 £270° 

and - 180 °<n ss -^< 0 o .*** As a result, the self-shadow interval is 

not necessarily symmetrical about i\ = 0 0 / as is the earth shadow 
interval C ± n__) . Consequently the self-shadow and earth shadow 

intervals may assume any one of four possible relationships, as 
shown in Figure (C-4) . 



Figure (C-4) - Possible Earth Shadow and Solar Panel Self-Shadow Regions (LV - Fixed Roll Mode) 


*See Footnote ** on page 45. 

**While the useful range on <j> for 6>0 is -90°j<<{><90 0 , it is 
only necessary in evaluating n gs ^ and n gg2 , to consider either 

-90 o <4><0 or 0<jj><90°. This follows from Eqs. (C-41) - (C-44) , 
which indicate that the extent of the self— shadow interval is 
the same, whether C>0 or C<0. The only difference is that for 
C<0 the self-shadow interval lies within the range 90°£ri£360 o 
whereas this range is 0°<_n£270° for C>0. Consequently only 
-90°_<(j)_<0 o is considered. 

***The convention adopted here is that n 0 is always a posi- 

ss ^ 

tive angle whereas n sg -^ is always a negative angle. 
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The terminal point (n ^.) at which the self-shadow 

region reduces to zero, occurs when n gs ^— ( ri ss 2 - ‘^^° ^ * ^is 

condition results when the radical in Eqs. (C-41) - (C-44) 
becomes zero or in view of Eq. (C-45) , when c$=ccipC4>. 

The self-shadow boundary is then specified by 


fl80° + tan 6 — cos 


n ss2 = < 


(C-46 ) 


no self-shadow 


6 > cos 1 (C 0 pCiji) 


and* 


r~90° - tan 1 (“4§~rr4|‘) £ — cos "^(ca^ciM 


L ssl 


no self-shadow 


6 > cos"*^ (ca c<j>) 
P 


( 047 ) 


and when S=cos (c a pC<f>), it follows from the preceding equations 
and Eqs. (039) and (040) that 


n sst n ssl 


= (n ss2 “360°) =-180° - tan 1 (-C/B) 
- -180° - tan 


•Y - . gM . t - \ # .|*| 

V*Ve 2 * - c 2 P,/ 


( 048 ) 


In Fig. (C-5) curves illustrating the behavior of n gg ^ 
and rt sg2 as a function of e and a p are shown for <j>= -45°. The 
terminator curve corresponding to n ^ is represented by the 
heavy dotted line. The earth shadow boundaries ( * ti ) corre- 

6S 

sponding to a 230 NM circular orbit are represented by the 
dashed curves. The various shaded regions in Fig. (C-5) denoted 
by A, B, C and D, correspond to the location of the self-shadow 


*This particular form for i gg ^ follows from tfie fact that 

,l s sr tanl ' 90 °' tan-lcc ' l ss i' /s ' l s si 1 1 ■ 
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Figure (C - 5) - Variation in Panel Self-Shadow Boundary with a p and /J (LV - Fixed Roll Mode: 0 = -45°). 
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a) 0=0° 
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b) 0= —90° 
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’fel > 

^ssl ^ ^ 


Figure (C - 6) - Variation in Panel Self-Shadow Boundary with dp and /? (LV - Fixed Roll Mode). 
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boundaries for the four cases illustrated in Fig. CC-4) . The 

relative size of these areas varies with. <{> and some vanish as 

<j>->0° or <j>+-90°. The curves for n. . and n. ~ and the corre- 

s s X ssz 

sponding shaded areas, when 4>=0° and <{>=-90 °, are shown in 
Figs. (C-6a) and (06b) respectively. When <{>=0°, the self- 
shadow interval is symmetrical about n=180°, as shown in 
Fig. (06a) where ri ?= -n ccl and r> . = -180° for all $ and a . 

SS^ SSI SSl p 

A 

Evaluation of E 

A 

The electrical energy (E) available from the solar 
array according to Eg. (5-16) is given by 


E 



r n 2 


C;)l LV dr| 


(049) 


where represents the array sunlight interval. In 

order to include the effect of distinct earth shadow and self- 
shadow intervals as shown in Fig. (04d) , it is convenient to 
modify this expression so that 


where 


E = 


2n 


es 


r 

r n 2 

CM 

\ 

cX LV d n - 

cX LV d n \ 

J 

J 

* 

^ n i 

t — i 
& 


(050) 


n 2 “ 


n l = 


ll ss2 

Tl ss2 <Tl es 

(Case : 

A) 


n es 

n es — ^ss2* 

(Case: 

B,C,D) 

(.05 la) 

^ssl 

^ ^ssl 1 <n es 

(Case: 

a,c,d) 


_t *es 

"es^KslI* 

(Case : 

B) 

(C-51b) 


*or if no self-shadowing occurs, i.e.. 


( 8>cos 


[ca c<f] ) . 

r 
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n 2 * 


n l “ 


^SSl 

n 360° <n. 

1 ss2 1 l es 

CCase: 

D) 

(c 

0 

otherwise** 

(Case: 

A ,B ,C) 

(n SS 2- 3 6°°) 

|„s 32 -36°°|<«e S 

(Case: 

D) 

(C 

0 

otherwise* 

(Case : 

A,B,C) 


After substitution for cA from Eq. (C-37) and integra- 
tion of Eg. (C-50) the expression for E becomes 


{[A(n 2 - 


n-^) + Bfsn^-sr^) + C (ct^-ch-l] 

- [A(n2~n 1 ') + B(sn 2 - sn 1 ) + C (c 




(C— 53) 

where A, 8 and C are defined by Eqs. (C-38) , (C-39) and (C-40) . 

A 

In Fig. (5-15) E vs. B is plotted as a function of the panel 
angle (a^) for three cases: <f>=0 , $=— 45° and <j>=-~90°. In each 

case the upper envelope boundary on the curves corresponds to 

A 

the maximum electrical energy (E ) available from the array 

max 2 

(for a specific <i>) with optimum variation of ot with £. Plots 

a ir 

of E and a , . . vs. $ obtained numerically for several 

max p (opt) 2 

values of <j> are shown in Figs. (5-16) and (5-17) respectively. 

These results indicate that making <J>=0° evidently yields the 

A 

largest E for all 8. Of theoretical interest perhaps, is 

IHdX 

the fact that in the region 0<B<23° other combinations of <t 

— ^ P 

and <j> also exist which yield the same as with 4>=0°, as 

Iuct.X 

will now be shown. 


*or if no self-shadowing occurs, i.e., (B>cos [ca ci)] ) 

ir 
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Ana lyt jcal Determination of ot^ and 

It is convenient at this point to introduce a change 
of variables for the terms involving a^ and 4> in the expression 

for cX^y in Eq.'(5-9). Specifically, let 


S(j>cctp = sj>_ sa 


(C-54) 


and 


sa = c* sa 

p x. - 


(C-55) 


where -90°< < 0° and 0°< a ^90°. Since 


2 . 2 


ca c<b — ■v/c^a - _ T _ „ 
p T V p T p 


~ *. n 2 t~ 

s 4> c = 'yc a^ — s ex 


2 2 2 

— *yc a + C j>_S a_ - s a_ = ca 


(C— 56) 


it follows that cX^y in Eq.(5-9) can be written as 


cX = ca c<j>s8 + sa cgcn - s^ca cgsri 
JjV P P P 

= ca_s8 + sacgc (j>+ri) 


CC— 57) 


The procedure here will be to express the electrical 

A 

energy E in terms of jx and j>_, next to determine the a and £ 

/v 

which maximize E and finally to infer from the relationships in 
Eqs. (C-54) , (C-55) and (C-56) , the corresponding optimum values 

A 

of a^ and > In evaluating E the panel self-shadow boundary is 

first expressed in terms of g, _a and from Eq. (C-57) after 
setting cX=0. This yields 


‘ c($+tj) = -tan8/tana_ 


8 <_ a 


(C-58) 
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so that 

. "SS2 = 5 " ± 

and 


where 


^ssl = "5 - i 


CC— 59) 


(C-60) 


5 


f 180° 
'.180° 


-1 


cos 


(tang/tanaj 


8 <a 


8>a 


For $=0 , the self-shadow boundaries are symmetrical 
about n=180°. Since £=0° implies from Eqs. {.C- 54) and (C— 55) 
that <f*=0° and a =a, the curves in Fig. (C-6a) for tin and n_o 

p *“ SS J. SS^ 

correspond directly to o ssl =“£ and n SS 2 =+ £ respectively. The 

effect of making -90° <£<0 is to shift these curves upwards, 
thereby creating four possible shadow boundary arrangements 
as described earlier in Fig. (C-4) . 

Expressions for the electrical energy (E) correspond- 
ing to each of the four shadow boundary arrangements in 
Fig. (.04) can be obtained after integrating Eq. (C— 50). Upon 
substituting, Eq. (C-57) for cX and appropriate values for the 

I I 

limits n^/ n^, and from Eqs. (C-51) , (C-52) and (C-60) 
the results are 


Case A: 


<- ri ss2 

[casg + sacgc (jf+n.) ] dn = ^ (casg • £ + sacgs?) 

n es 

a ssl 


(C— 62) 
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Case B: 


E = 


2n 


es 




[c£Sg + sacgc U+nJ Jdn = (casg»rL es + S£Cgsn es c<}>) 


‘es 


- n 


es 


Case C: 


(C— 63) 


n 


es 


E = 


2n 


tcasg + sacgc((j)+n) ]dn 


es 


n 


ssl 


= 2r^ {“ ! 

'es V. 


casg (n_ +<f>+£) + sacg Is (n. _+<(>) + s£;] 

co 65 — 


} 


(C-64) 


Case D: 


E = 


/- n es 
f es 

^ssl 

cX LV dn - 

oA LV an 

1 J _ , J 

\ 


es 


ss2 


1 ft 

n es l " 


C£s8*ri + sacBsri 0 c<j>] - [casg (2 tt-?) 

CO 65 "** 


- sacBsS] 


■) 


(C-65) 


The optimum a and £ can be obtained in each case by 
setting 9E/3o=0 and 3E/3£=0 respectively. In evaluating the 

results however, only those solutions for a and ^ . which 

opt Opt 
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yield shadow boundaries consistent with the respective cases 
CA, B, C or D) are acceptable.* 


Case B: 


and 


It follows from Eq. (C-63) that & O p t = 0° 


3 E X 

jg - — (-sasB-n es + c«c6sn es ) - 0 

6S 


so that 


tan ^ P t 


STl es //r| es 

tang 


and 


(C-66) 


(C-67) 


ca 


p Copt) 


tang 


>/csn es /n es )^ + tan 2 6 


(C-68) 


Consequently E max in case B is given by 


max 


" ^ {w Se %*t[ l + ( SfL tane eS ) tan 2opt]} 


- ( £e/c ^pt) = 4 


c 2 gs 2 rL es + s 2 B ♦ n^ s CC-69) 


*With regard to case D however, it follows from Eqs. (C-63) 
and (C-65) that E^ d j<E^ b ^ so that no further consideration of 

this case is actually necessary. 
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Hence, 


It follows from Eq. (062) that E is independent of £. 


9E _ 

Sot " T) 


es 


3 £ 

-sasg*C + cacgsg + (casg + sacgc£)-r-2. 
— — — —.9a 


= 0 (070) 


so that 


sa sg*£ = ca eg s £ 


(071) 


since 


c£ = tang/tana 


(072) 


from Eq. (061) . Since ( is a function of £, it is necessary to 
eliminate a via Eqs. (071) and (072) , solve for £=£ 0 p t as a 

function of g and then determine ^ from Eq. (071) where 


tana . = 3S °P t/ ^°P t 
—opt tang 


(073) 


Eliminating £ from Eqs. (071) and (072) yields the 
following transcendental equation 


2 £ = - 


s2? 

2 

tan g 


(074) 


As illustrated m Fig. (07a) two solutions for £ exist for 
sufficiently small g. The limiting case corresponding to g=g c 




Figure (C - 7) - Relationship between Panel Self-Shadow Parameter { |) and (3- 
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occurs when c2£= -tan g which means, in view of Eg. CC-74) , that 


2£ c = tan2£ c 


(C— 75} 


Solution of Eq. (C-75) yields £ c =128.7°, which implies that 

8 c =tan ^/-c^2£ c }=25. 0°. The solutions for £ = 5 Q p t obtained from 
Eq. (C-74) and a curve of the earth shadow boundary (n ) are 

6S 

plotted in Fig. (C- 7b) as a function of 6. The intersection of 
the £ and n curves represents the transition point from 

Opt c& 

case A to case B. In view of Fig. (C-7b) and Eq. (C-59) , case A 
applies provided 


or 


min 


n ss2 ^opt — n es 
5 ( - n es + W £ ♦ £ °° 


(C-76) 


(C— 77) 


For 8=0° it follows that <6 . =C-n. +90°) , since £ .=90° whereas 

—mm es opt 

tnirT 0 at 6=B es' since «opt =, W 


The optimum a, for case A is given by Eq. (C-73) and is 
similar in form to the result in Eq. (C-67) for case B. The 

A 

corresponding electrical energy (E ) obtained from Eq. (C-62) 
with a=a 0 p t is given by 


max n 


1 7 ~ 2 n .~2 


es 


c es £ opt + s 6 * 5 opt 


CC— 78) 


which is identical with Eq. (C-69) when £ =n 

opt es 
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Case C : 

A 

The procedure is similar. Setting 3E/3£=0 results in 


^opt 


C-a es+ ?) 


CC-79) 


Substituting this result into Eq. CC-64) yields an expression 

A 

for E which is identical to E for case A in Eq. (C-62) . The 

result obtained from SE/Da^O is also identical to Eq. CC-71) so 
that case C is really included in case A insofar as determining 
the optimum ct and is concerned. 

In summary then, it follows that 


°±^ es : 


&es- B - 90 ° 


E = -i- V c2 ps 2 £ + s 2 g*g 2 . 
max xi opt opt 

65 


tana - S *°P t ^ qpt 
—opt tang 


‘-■W^opt 1 iii°' 


c = -i- *>/c 2 gs 2 n. + s 2 g*rt^ 

max it „ v y L es y L es 

es 


tana , 
-opt 


srL es //,1 es 

tang 


iopt 


CC— 80) 


CC-81) 
(C— 82) 


CC— 83) 


CC— 84) 


0 ° 


CC— 851 
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The optimum spacecraft roll angle C<J>) and optimum panel angle 
(a^) corresponding to and are than obtained from 

Eqs. (C-54) , CC-55) and (C-56) , which can be expressed as 


tan<{> = s<j) tana 


(C-86) 


and 


sa 

P 


= Cj>_ so 


CC— 87) 


Curves for a p (opt) an< ^ ^opt obtained from Eqs.(C-86) and (C-87) 

for iopt =0 ° and ^opt =(_n es +C opt ) with ^opt given by Egs - (C “ 81) 
and (C— 83) are shown in Fig. (C-8) . In the region 0<&<8 =22.9° 

@ S 

these curves represent the bounds on the other possible combi- 

A 

nations of a , , . and <f> . which yield the same E . For 

p (opt) opt * max 

£=0° it follows from Eqs. (C-86) , (C-87) and (C-84) that 4>=0 ° 

and a , =a . . 
p(opt) ~opt 


W (b) 



Figure (C-8) - Optimum a p and <p vs j6 (LV - Fixed Roll Mode) . 
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C.5 LH-Optimum Roll Mode 


Panel self-shadow boundaries in LH inodes are closely 

related to those in corresponding LV inodes. The expression for 

cA T „ in Eq. (5-41) with defined by Eq. (5-43) is 
Lii 


CA 


LH 


= Cl 


[- sa I 


an + co 


■v/tan^B+c^ 


:] 


(C-88) 


If n is replaced by (n-90°) it follows that Eq. (C-88) is iden- 
tical to cA t , 7 in Eq. (C-3) . Hence self-shadow boundaries have the 

same general character as in the LV-Optimum Roll mode, but are 
symmetrical about 90° rather than 180°. In view of Eq. (C-6) 
the self-shadow entry and exit points (n ss2 • n. ss i) are 9^ ven 


"ss2 - 180 °- T 'ssl 


^90° - cos (co /c8) 6 £ a 


(^no self-shadow 


> a 

P 

(C-89) 


A plot of (n ~,n ,) vs. 8 is shown in Fig. (C-9) where the 

SSb S S i. 

dashed curve represents the earth shadow boundary (n eg ) • 



Legend: 
Case C 


T?SSl >% 


es 


Case D □ ^ssl 


Figure (C-9) - Variation in Panel Self-Shadow Boundary with a p and 0(LH - Optimum Roll Mode). 
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Self-shadow and earth shadow regions are distinct, 
whenever n. n <;rL (.Case Dl and overlap partially whenever 

S.S 1 6S A 

n. i >n (Case C) as shown in Fig. (5—20) . To evaluate E it is 

SSI 6S 

necessary to modify Eq. (5-16) such that 



(C-90 ) 


where 



(C-91a) 



(C-91b) 



(Case C) 


(Case D) 


(C-91c) 


Rearranging cX similar to the procedure in Eq. (C-8) and sub- 
LH 


sequent integration of Eq. (C-90) yields* 


(Case C) : 



(C-93a) 


(Case D) : 



(C-93b) 


where 


k = c6 


(C— 9 4 ) 


*See Footnote *, p. C-4. 
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Formal analytical evaluation of E with respect to a /as in 

max p 

the LV modes, is not amenable to closed form solution, since E 
is a more com] 

Nevertheless, E_ 


is a more complex function of a in both case C and case D. 

P 

max and the corresponding a p£ 0 p t ) were obtained 

numerically as shown in Figs. (5-29) and (5-30). Since 0t p( O pt) = ®° 

for all 6, it follows from Fig. (C-9) that case D applies and 
n gS 2 = 90 o . Hence, Eq. (C-93b) yields 


E = E(k, 
max 


n es^ //n es 


(C-95) 


C.6 LH-Cyclic Roll Mode 

In the cyclic roll case, the expression for cA TtJ in 

Ln 

Eq. (5-41) with <f> defined by Eq. (5-44) is 


(CX LH»1 " 

/ v 2 

-sns(a + 6) + ca cgc n 
P P 

-180°_<n<0 o 

(C- 

■95) 

(CX LH*2 = 

, x 2 

-sns a -3) + ca c£c n 

p p 

0° <^n^_180° 

(C- 

-96) 


If n is replaced by (n-90°) it follows that Eqs.(C-95) and 
(C-96) are identical to cA^y in Eqs.(C-26) and (C-27) . Hence 

self-shadow boundaries have the same character as in the LV- 
Cyclic Roll mode, but are symmetrical about 90° rather than 
180°. In view of Eq. (C-29) the self-shadow entry and exit 
points < I ' SS 2 ' ’’ssl* are 9iven by 





— 

1 

f 90 °-cos -1 

-s(a p -8) + ■ 

■\js 2 (a p -6) +4c 2 a p c 2 6 

180 °-n ,= 

ssl 

1 


2ca c6 

P -1 

1 

i.no self-shadow 





6><, p 


(C-97) 
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A plot of Cn ss 2 • ^ssl^ vs * ^ sflown i- n Fig. CC-10) where the 
dashed curve represents the earth shadow boundary (n ) . 

6 S 


(deg.) 



0 15 30 45 60 75 90 

P( deg) 


Legend: 



Case C 

■ 

^ssl ^^es 

Case D 

□ 

^ssl ^^es 


Figure (C - 10) • Variation in Panel Self-Shadow Boundary with a p and /J (LH - Cyclic Roll Mode). 


Self-shadow and earth shadow regions are distinct, whenever 
ri „ (Case D) and overlap partially whenever n _-,>n 

SSI 6S SSI 6 S 

(Case C) as shown in Fig. (5-20). The modification of Eq. (5-16) 

A 

to evaluate E is 


E = 


2 n 


es 


f 


r 2 

CM 


^ CX LH^l dn + 

^ C ^ LH ^2 

2^ n 


j 

n * 

' J 


(C-98) 
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where 


= -'ll, - n. ; 


es 


(C-99a) 


n l = n ss2 


CC-99b) 


I 



r n Case C 

es 

.. n ssl Case D 


(C-99c) 


Integration of Eq. (C-98) after substitution for ( C X LH ) ' j_ and 
(ex 2 yields 


(Case C) : 

e = 2rieg ( 2c V 8 "* 3 (ctp+6) c n es + s ^p “^^ Cr > SS 2 


» ca eg r .'x 

+ L 2(t W n ss2 ) + (s2n es +s2n ss2 ) J) CC " 100a) 

(Case D) ! 

E = fv 6 ( 1 - cn es> + s( V 3)cT W 

ca eft f > 

+ -4 L 2C W n SS 2-" /2) + (s2ri es +s2T1 S s2 ) ]j (C-100b) 


Again formal analytical evaluation of E m ^ y with respect to 
is not amenable to closed form solution. Numerical results 
for E max and a p( Q p t j i n Figs. (5-29) and (5-30) indicate that 
except for =10°, E__„ is obtained for a. ,• ..=0°. Prom 

Pig. (5-22) it is clear that for a^O 0 , E is also close to the 
optimum for small ft . For a p~0° case D applies and 9 gS 2~90° as 



bellcomm, Inc. 

Appendix C 


C-34 


indicated by Fig. (C-10) . Hence, Eq. CC-100b) yields 

E max " [sBCl-ci es ) + =6 C2a es +s2a es )/4]/ %s 6 > B cyc - 10“ 

(C-101) 

For $<3 =10°, this result is also a reasonable approximation 

cy c 

to the optimum. 

C.7 LH-Fixed Roll Mode 

In the fixed roll case <j> is constant m Eq. (5-41) , 
which can be written as 


cA t „ = A - Bsn - Ccn 

hn. 


CC-10 2) 


where A, B and C are as defined in Eqs. (C-38) , CC-39) and 

(C-40) . If n is replaced by (n-90°) , this expression is iden- 
tical to cA^ in Eq. (C-37) . Hence for a given 4> , a and 8r 

-uV P 


the self-shadow boundary points, n ssl and n ss2 are merely rotated 

by 90°, clockwise in Fig. (C-4) . This yields the four possible 
shadow arrangements, cases A-D, shown previously in Fig. C5-23) . 

In view of Eqs. (C-46) , (C-47) and (C-4 8) , the corresponding 

expressions for n o1 , n 0 and n . , the self-shadow terminal 

SS J. SS4 SSt 

point, are given by* 



\ no self-shadow 


6 

3 


< 


> 


-1 

cos 

-1 

cos 


(ca c<j>) 

ir 

(ca cd> ) 
P 


(C-103) 



v no self-shadow 


8 

8 


< 


> 


-1 

cos 

-1 


cos 


(C 0 tpC<}>) 
(c« Cf) 

ir 


(C-104) 


*See Footnote ***, p. c-15. 
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and* 

n sst = 9Q ° ~ tan" 1 (c g s §/" m \[c 2 <f> - c 2 8 ) |g| > | | CC-105) 


In Fig. (C- li) curves of n , n.., and n and the four possi- 

ble regions corresponding to self-shadow cases A-D, are shown 
for <j>= -45°. The relative shape and size of these regions 
differ from the LV-Fixed Roll mode," but also vary with <f> with 
some vanishing as <j>->-0 0 or 90°. The curves for n ssl and n ss 2 

with <j,=0° and ^=-90° are shown in Fig. (C-12) . 

/s 

Evaluation of E based on modifying Eq. (5-16) , as in 
Eq. (C-90) , with Eq. (C-102) substituted for cA^ yields 

~ i ( 

E = jz — < [Atru-n-j) + 8(sn 2 _si i 1 ) - 0 ( 0 : 12 - 012 )] 
n es V. 

ft It 

- [A (112-02) + B(sn 2 “ s n 2 ) 


- C (cn 


2' cri i 1 } 


(C-106) 


*While (o gst ) L Y in E q. (C-48) was defined as a negative 

angle, a positive convention was used in Eq. (C-105) since 
90°_< (n gst )f,180° . This represents 


* n sst^LH ^ n sst^LV 


90 0 ] + 360° 
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(deg.) 
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Legend: 



Case A 

■ 

^ssl ' > ^es 
*?ss2 

Case B 

□ 

^ssl ^ -T ies 
^ss2 ^ ^es 

Case C 

■ 

^ssl <- ^es 
r ks2 <r kt 

Case D 

□ 

T? ss2 <T ?es 
r?ss1 + 360 <7j ( 


Figure (C - 11) 


Variation in Panel Self-Shadow Boundary with d p and ft (LH - Fixed Roll Mode: 0 = -45°). 
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’fei 

(deg.) 


(a) <t> = 0° 
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Legend: 
Case C 


1 \s2 <7 ?« 


Case D I 1 iJ ss j +3 60°<r) eJ 


Legend: 


Case A 

nssi>-nes 
*1* 2 <r >es 

CaseC 

nssi<-nBs 
1 ?ss2 ^es 


Figure (C - 12) Variation in Panel Self-Shadow Boundary with Op and /3 (LH - Fixed Roll Mode: 0 = 0°, 0= -90°). 
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where n 


are defined by 




f ^332 

n ss2 < "es 

CCase: 

A f C/D) 

n 2 = ( 

I 

l V 

"as i "s=2* 

(Case: 

B) 


' "ssl 

l n ssll ' "es 

(Case: 

A f D) 

11 1 =< 

, r 'es 

l"ssJ - "es* 

(Case: 

B,C) 

f 

f(n S sl +360 °) 

Cn ssl+ 360«) < n es 

(Case: 

D) 

n 2 = * 

) 

0 

otherwise* 

(Case: 

A/B,C) 

t 

' n ss2 

<"ssl +36 °°> < "es 

(Case: 

D) 

n l = < 

0 

otherwise* 

(Case: 

A/B,C) 


(0107a) 


(0107b) 


(C-108a) 


(C-108b) 


Determination of E. 


max 7 a p(opt) 


and 4> 




The procedure and results obtained in Section C.4 for 
the LV-Fixed Roll mode can be used with slight modification to 


o 


obtain E max and the corresponding optimum a p and 

to Eqs . (C-54) and (055) the change of variables 
the LH case is 


(f> . Analogous 
introduced in 


sd>ca = — c<f> sa 
P x “ 

(C-109) 

sa = s<f> sa 

p - - 

(C-110) 

where again 0<ct<90°/ but 0<£<90°. Since 

co cd> = ca 
P - 

Cc-111) 


(f?>cos 1 [c apC( } l ] ) . 


*0r if no self-shadowing occurs/ i.e., 
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it follows that cX LH in Eg. (5-41) can be written as 


c\ TTJ — ca sg + sa eg c(d»+n) 

Juii — — — 


(C-112) 


which is identical to Eq. (C-57) . Consequently self-shadow bound- 
aries can also be expressed as in Eqs. (C-58)- (C-61) . The result- 
ing evaluation of E for cases A and B is identical to Eqs. (C-62) 


and (C-63) . 
ent of <f> . 


For case B, an <3 for case A,. 


E is independ- 
max 


Determination of a p ( opt ) leads to Eqs. (C-67) and (C-73) 

where ? ^ is the solution of Eq. (C-74) shown in Fig. (C-7) . The 

transition point from case A to case B corresponds to the inter- 
section of the C and n curves at g=g =22.9° for a 230 NM 

opt es es 

orbit. In view of Fig. (C- 7b) and Eqs.(C-60) and (C-107b) , case 
A applies provided 


or 


^ssl ^opt — > ^es 


0 < A < n — £ . = 4 

x es opt -Hnax 


(C-113) 

(C-114) 


E , 
max 


2opt and ±opt are 


In summary then , the optimum results , 

identical to those in Eqs. (C-80 ) - (C-85) for the LV-Fixed Foil 
mode except that Eq. (C-114) replaces Eq. (C-82) for the condition 
on 4 . 


The optimum spacecraft roll angle (<J>) and panel angle 
(cip) corresponding to 4^ ^ and ^ are obtained via Eqs. (C-l-09)- 

(C-lll) which can be expressed as 


tan<j> = -c4 tan£ (C-115) 

so “ s_4 sa (C-116) 

Curves for a p( 0 p t ) and ^ 0 pt obtained from Eqs. (C-115) and (C-116) 

for iopt =0 ° and lopt =n es _5 opt with ^opt given by Eqs - (C " 81) and 
(C-83) are shown in Fig. (C-13) . In the region 0<g<g =23° these 

C S 

curves represent the bounds on the other possible combinations 
of °p ( 0 pt) and ^opt yield the same E max - For £=0° it 

follows from Eqs. (C-115) and (C-116) that a , ^.=0° and 4 , = 

p (opt) opt 

-a . This indicates that -4 .in the LH-Fixed Roll mode is 
—opt: opt 

equal to a^ ^ ^ in the LV-Fixed Roll mode. 
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Pointing Angles for Sun-Oriented Solar Arrays and Experiments 
on Spacecraft in LV, LH and POP Modes 


In this appendix pointing angles for sun-oriented solar 
arrays and experiments on spacecraft in LV, LH and POP modes are 
evaluated from appropriate coordinate transformations. 

D. 1 Solar Array Pointing 


The coordinates (x ,y ,z), Cx ,y ,z ) and {.x ,y ,z ) 

p p p i\ T[ J\ boo 

shown in Fig. CD-I) are defined as the solar array, orbit and 

sun referenced coordinate systems. From the transformation 

relating tx .y.z ) and Cx.y.z ) for LV, LH and POP space- 
P P P s s s 

craft attitude modes it is possible to determine: 


1) 


2 ) 


the pointing angles (cip,<}>) for aligning 

the outward normal to the panels (z^) with 

the sun line Cz ) and 
s 

the incidence angle (X) between the sun line 

Cz } and the panel outward normal (z ) for 
s p 

arbitrary oip,<{>,ri and g. 







BELLCOMM, INC. 

Appendix D 


D-2 


D.1.1 Transformations 


The transformation relating (x ,y , z ) with. 

n ■‘n n. 

( W z s> is 


% 


= T Y T 


X 


n -6 


-B ^ 


(D-l) 


The transformation relating (x ,y , z ) with (x ,y ,z ), 

P P P po ,j: po po 

which represents the panel coordinates with a =$=0 , is given 

P 

by 


_ T Y T x 
-P a p <t> *P° 


(D-2) 


The transformation relating (x .y . z ) with (x ,y ,z ) is 

po po po n ,J n n 


u po = B ± y n {i = LV , LH , POP ) 


(D-3) 


where is to be specified for each attitude mode. From 
Fig. (D-2) it follows that* 



Figure (D-2) - Orientation of Solar Array Coordinates in LV, LH and POP Attitude Modes with a p = 0= 0°. 


*See also Figs. (5-7) , (5-18)' and (5-27). 
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R = T 1 ^ m Y 
*LV 90 -90 


0 

-1 


0 

0 


0 -1 


1 

0 

0 


CD-4) 


and 


R = m x 

LH 90- 


P — T 1 

POP 90 


1 

0 

0 


0 

-1 

0 


0 

0 

-1 


1 

0 

0 


0 

1 

0 


0 

0 

1 


(D-5) 


(D-6 ) 


Cons equent ly 


y = A . y 
-=4? i — s 


(D-7) 


where 


A. = T y Tf B T y T * (i = LV , LH , POP ) 
l a (J> 1 n - e 


CD- 8} 


D.1.2 Pointing Angles (a p/ <j>) 


With the panels pointed at the sun the solar vector 


is 




= A . S_ = A. 
l— s x 



(D— 9 ) 



Bellcomm, Inc. d-4 

Appendix D 


or in view of Eqs.CD-7)and (D-8) 



(LV Mode ) 

sa^ = cncg 
tan^i = -sn/tanp 

(LH Mode ) 

s a = -sncg 
P 

tan^ = ~cn/tang 

CPOP Mode ) 

sa = -sB 

ir 

tan<}> = -tann 


CD 

CD 

CD 

(D 

(D 

CD 


D.1.2 Sun Incidence Angle (A) 

With the panels not pointed at the sun the solar 
vector expressed in ( x p'Yp' z p) is 



The incidence angle CA^) is then specified by 


11 ) 
- 12 ) 

-13) 

-14) 

-15) 

-16) 


CD-17) 
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D-5 


so that 



CD- 19) 


(D- 20) 


CD-21) 


D. 2 Experiment Pointing 

The results in Section D.l for solar arrays are also 
applicable to determining pointing angles for sun-oriented experi- 
ments on spacecraft in LV, LH and POP modes. The experiment 
degrees-of-f reedom are assumed to be a rotation (<f> e ) about the 

spacecraft roll axis and a rotation (a g ) about an axis normal to 

the roll axis.* With “ e = 4* e = 0° and no spacecraft roll (<j>=0°) 

the experiment pointing axis is assumed parallel to the spacecraft 
roll axis, z q/ shown in Fig. (D-l) . Hence, a g and <j> e are equiva- 
lent to a and <j> in Eqs. (D-ll) through (D-16).** To account for 
P 

cases where (j> ^ 0°, Eqs. (D-12) , {D-14) , and {D-16) are merely 

modified such that C^ e -4>) • 


*The rotation <f> could be realized partially with internal 
e 4 

gimbaling or exactly by a spacecraft configuration comprised of 
modules stacked along the roll axis, each with an independent roll 
capability. 

**As a practical matter experiment occultation by the spacecraft 
may occur when the pointing axis is oriented nearly parallel to the 
roll axis unless the experiment module is located at one end of the 
spacecraft to permit axial pointing. 
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Sunlight Interval in Elliptical Orbits 


The geometry' for determining the sunlight interval in 
elliptical orbits is similar to that used for circular orbits 
as shown in Fig. (E-l) . Here however, the orbit altitudes at P^ 

and P 2 differ as given in Fig. (E-lb) where ri^ and fj 2 represent 

the angular displacement of P 1 and P 2 from perigee and n s is 

the location of midnight relative to perigee. The relationships 

13 

between (o^,a 2 ) and g are given by 


C0^ = cp 1 /cg 

(E-l) 

C 0 2 = cp 2 /cg 

(E-2) 

sp-^ = Sff [1+ec (ti s -0 1 ) ] 

(E-3) 

s p 2 = S0 [1+ec (n s +0 1 ) ] 

(E-4) 

sp — R/a(l-e^) 

(E-5) 


where R is the earth's radius, a the orbit semi -major axis and 
e the orbital eccentricity. In view of Fig. (E-2) , n g is given 

by 


n = ifi+180 °-ft • (E-6) 

S a. 

where S2 is the location of perigee relative to the ascending 

cl 

node as defined in Eq. (2-11) . 

For a given g and n s , solution of Eqs. (E-l) through 
(E-4) yields a-^ and 0 2 « The time in shadow obtained from Kepler's 



BELLCOMM, INC. 
Appendix E 


E-2 


equation in terms of the eccentric anomalies E^ and E 2 is 


II 

d 

U) 

Eh 

Ji I03 2 -esE 2 ! - ■ 

(E— 7) 

\ 

SE i = 

/l-e 2 sn 1 

CE-8) 

1+ecfi^ 

SE 2 - 

L 2 - 

Jl-e sri 2 

CE-9) 

l+ecn 2 

*1 = 

n s‘°l 

CE-10) 

~ n 2 = 

V°2 

Ce— 11) 


and T is the orbital period. The sunlight fraction of the 
orbit is then 

T SL - 1 - T sh /T - [< e 2- E 1> - stsE 2- SE i n CE - 12) 


For circular orbits, e=0 so that a=R+H, p^=p 2 =0 an< ^ a l =0 2* 
Thus , 


ir- a 


sL 


= 1 - 


2 it (2a l ) 


1 - 


es 


(E-13) 
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where r is one-half the orbit sunlight angle in radians, 

65 


and 


Cct-l = c<j/c8 * (E— 14) 

so = R/CR+H) CE-15) 


*While umbra and penumbra aspects of the earth-sun geome- 
try have not been included here, this slight effect can be 
included easily. The shadow interval boundaries are based on 
visibility of the full solar disk for the penumbra effect and 
the first (or last) point of visibility of the edge of the 
solar disk for the umbra effect. In Eq. (E-14) , for example, 
cct becomes c(cr+a ) for the penumbra effect and c (o-a ) for the 

umbra effect where- a s = 0 . 26 ° is half the angle subtended by the 

sun at earth distances (1AU) .' The same modification applies to 
cp^ and cp£ in Eqs. (E-l) and \(E- 2 ) . 
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List of 


For convenience of the reader this list contains all 
important symbols, a brief definition and the page location 
where first used. Minor symbols appearing in only one section 
are omitted. 


Greek 

Symbols 


Definition 


First Used 
on Page 


sun line declination relative to equatorial 
plane 

launch azimuth relative to geographical 
north 

solar array articulation angle 


auxiliary variable used with £ to define a 
in LV and LH-Fixed Roll modes ^ 


C-21 


p (opt) 


optimum array articulation angle 


a . 
—opt 


optimum value of a 

minimum angle between sun line and orbital 
plane (sun line declination relative to 
orbital plane) 

angle denoting 8 envelope boundaries 

constant angle about which 8 varies in sun- 
synchronous orbits 

maximum possible | 8 | for particular orbit 
inclination 

dg/dt 


C-27 


'max, max 


maximum of | 8 | for particular orbit inclina- 
tion and h 

rm 

maximum of |3[ max With respect to orbit 
inclination 


angle denoting earth- sun line location in 
ecliptic plane relative to autumnal equinox 
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Greek First Used 


Symbols 

Definition 

on Page 

y 

! o 

Y at time, t=0 

13 

Y 

r when IsHsInax 

18 

• 

Y 

mean rotation rate of earth-sun line about 
the sun 

7 

e 

eccentricity of satellite orbit about earth 

10 

e 

s 

eccentricity of earth’s orbit about the sun 

7 

n 

angle denoting spacecraft location in orbit 
relative to orbital noon 

41 

A 

n 

angle denoting spacecraft location in orbit 
relative to ascending node 

35 

*1 

n when array exits from shadow region 
(earth or self-shadow, circular orbits) 

42 

^2 

when array enters shadow region (earth 
or self-shadow, circular orbits) 

42 

^es 

n when spacecraft enters earth shadow; 
one-half orbit sunlight interval (circular 
orbits) 

37 

n eso 

n at 8=0° 

es 

44 

n ssl 

n when array exits self-shadow region 

42 

n ss2 

t| when array enters self-shadow region 

42 

fl l 

spacecraft true anomaly when entering earth 
shadow 

E-l 

fi 2 

spacecraft true anomaly when exiting earth 
shadow 

E-l 

n s 

location of orbital midnight relative to 
perigee 

E-l 

X 

sun line incidence angle on solar array (in 
LV, LH and POP modes: ^POP^ 

41 

x_ 

maximum projected area ratio between ends 
and lateral surface y ofL_gy linder/hemi spheri- 
cal-end spacecraft 

77 
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Definition 

earth's gravitational constant 

general unit vector expressed in coordinate 
system (i) with components e.g. 

li and ji n for unit vectors in spacecraft 

and nodal coordinate systems 

spacecraft roll angle 

<J> in LV-Continuous Roll Mode 

<j> in LV and LH-Cyclic Roll Modes 

optimum value of <{> in Fixed Roll modes ($=0) 

auxiliary variable used with £ to define 
in LV and LH-Fixed Roll Modes 

optimum value of £ 


First Used 
on Page 


C-21 


C-27 


angular displacement of spacecraft principal 
and geometric axes about roll axis 

dij>/dt 

angle between orbital noon and ascending 
node (sun line right ascension relative to 
ascending node) 

difi/dt 

average of iji over time 
maximum of i over time 

orbital angular velocity in circular orbits 
(constant) 

magnitude of bias gravity gradient torque - 
impulse per orbit (bias momentum) 

angle between perigee and orbit ascending 
node 

angle between orbit ascending node and 
autumnal equinox (measured in equatorial 
plane) 
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Greek 

Symbols 


po 

£2 


SI 


r 


& 


r 


Si 


L 


SI 


N 


rm 


Definition 


First Used 
on Page 


£2 at time, t=0 

d 

£2 at time, t=0 
r ' 

£2 r when Nl = |el max 
(£2 rQ -Y o ) , a constant 

launch site longitude relative to orbit 
ascending node (measured in equatorial 
plane) 


9 

9 

18 

29 

14 


longitude of equatorial noon meridian plane 
relative to autumnal equinox (measured in 


equatorial plane) 14 

rotation rate of line-of-apsides 9 

rotation rate of orbit ascending node 9 

£2„ when i=0° 10 

r 


English 

Symbols 


a 

a 

e 

f (3) 
f av 


<f > 
av 

<<f >> 

av 


orbit semi-major axis 
earth albedo constant 

fixed angle (23.45°) between equatorial and 
ecliptic planes 

functional relationship between a system 
performance factor and 8 

mission average of f(3) 

mean of f relative to launch time of day 

clV 

(t L ) and launch date (D^) 

minimum of f „ relative to t T and D 

av L L 

maximum of f relative to t T and D T 

av L L 


10 

73 

8 

82 

82 

83 

82 

82 


i 


orbit inclination 


8 
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English 

Symbols 

Ss 

t 



First Used 


Definition on Page 

i for sun- synchronous orbits 20 

launch azimuth 8 

- current time measured from a prior event 
(e.g. launch, orbital noon, autumnal equinox, 
etc.) 7 

launch time of day at launch site (U.T.) 13 

orbit time corresponding to 42 

orbit time corresponding to 42 


A maximum projected area of a geometrical 



surface 

73 

D 

diameter of a cylindrical spacecraft 

77 

d l 

launch date at launch site (calendar days 
since 00^:00 U.T. of September 23) 

13 

E 

electrical energy per orbit generated by a 
solar array 

42 

E 

s 

E with the array in a sun-oriented mode 

43 

E 

so 

E when p=0° 
s 

44 

A 

E 

A 

E normalized to E (E=E/E ) 

s s 

44 

A 

E o 

A 

E normalized to E (E =E/E ) 

so v o 7 so' 

44 

E max 

E maximized with respect to 

48 

E 

o max 

A 

E Q maximized with respect to a 

A A A P 

(E ~E E ) 

o max s max' 

68 

A 

E 

s 

A 

E normalized to E (E=E/E) 

s so s s so 

44 

F 

DS 

view factor for direct solar radiation to a 
surface 

73 

F 

ET 

view factor for earth-emitted radiation to a 
surface 

73 
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English 

Symbols Definition 


First Used 
on Page 


F 


RS 


view factor for reflected solar radiation 

to a surface 7'3 


F DS + aF RS 


74 


H 


circular orbit altitude 


10 


L 


length of a cylindrical spacecraft 


77 


F i 


unit vector specifying pointing axis to be 
aligned with the solar vector* 11 


P 



instantaneous power generated by a solar 

array 41 

P with the array in a sun- oriented mode 41 

earth's radius (3443.9 NM) 10 


S. 
— i 


solar vector, a unit vector* along the earth- 

sun line directed toward the sun 6 


S 


T 



T 


SL 



rnY 

V T ? 


solar constant (443 BTu/hr/ft 2 ) 12 73 

orbital period (2ir/o) o in circular orbits) 35 

gravity gradient torque (vector*) 35 

mission duration 82 

fraction of an orbit in sunlight 39 

elementary coordinate transformation matrices 5 

(see definitions in Section 2.1) 


V7 


n 


W. 


DS 


W. 


ET 


direct solar component of incident thermal 
energy per orbit with A normal to the solar 
vector 

direct solar component of incident thermal 
energy per orbit on a surface (normalized to 


earth-emitted component of incident thermal 
energy per orbit on a surface (normalized to 


75 


74 


74 


*Subscript i denotes coordinate system in which the vector 
is expressed. 
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English 

Symbols Definition 


W 


RS 


W 


s 


reflected solar component of incident 
thermal energy per orbit on a surface 
(normalized to W ) 



+ W 


RS 


First Used 
on Page 


74 

77 


Miscellaneous 
Symbols 


magnitude of a scalar or vector 


3 


t 


matrix transpose 


5 
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